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CHAPTER  1 
INTRODUCTION 


1.1.  A  Brief  Overview  of  Descriptor  Variable  and  Singular 
Perturbation  Theory 

Descriptor  variable  theory  depends  heavily  on  the  theory  of  regular 
matrix  pencils  which  was  first  explored  by  Weierstrass  in  the  nineteenth 
century  [1].  A  standard  modem  reference  to  the  theory  is  Gantmacher  [2]. 

A  regular  matrix  pencil  is  a  matrix  polynomial  Es-A  where  E  and  A  are  square 
matrices  of  the  same  dimensions  and 


det(Es-A)  $  0. 


(1.1) 


In  (1.1)  the  determinant  is  formed  in  the  obvious  way  by  taking  the  deter¬ 
minant  of  the  corresponding  matrix  of  scalar  polynomials. 

More  recently,  it  was  observed  by  Rosenbrock  [3], [4]  that  the 
linear  system 


Ex  =  Ax  +  Bu 


(1.2) 


is  strongly  related  to  the  pencil  Es-A  since  Laplace  transformation  of  (1.2) 
yields 

(Es-A)x  =  Bu  +  Ex(0) .  (1.3) 

x  and  u  are  the  Laplace  transforms  of  x  and  u.  If  E  is  singular,  (1.2)  is 
called  a  descriptor  variable  system. 

In  [4]  Rosenbrock  introduced  his  decomposition  of  (1.3)  into  static 
and  dynamic  parts  along  with  the  theory  of  infinite  decoupling  zeros.  The 
decomposition  uses  the  canonical  form  of  a  regular  pencil  to  decompose  the 
system  into  two  parts,  one  whose  eigenvalues  are  the  same  as  those  of  (1.2) 
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and  one  with  no  dynamics  in  the  usual  sense.  The  theory  of  infinite 
decoupling  zeros  is  an  algebraic  characterization  of  controllability  and 
observability  for  the  static  subsystem.  The  theory  utilizes  properties  of 
the  pencil  As-E.  Continuing  in  the  same  direction,  Verghese  et  al.  [5]- [7] 
carried  the  work  of  Rosenbrock  further  with  some  modifications. 

Proceeding  in  a  somewhat  different  direction,  Luenberger  (who 
originally  coined  the  phrase  "descriptor  variable")  concerned  himself  wii 
time  varying,  discrete-time  descriptor  systems.  In  [8]  and  [9],  as  well 
having  stated  basic  results,  he  gave  many  real-world  examples  as  a  justifi¬ 
cation  of  the  development.  [10]  contains  applications  of  the  theory  to 
the  LQ  regulator  problem  and  to  large-scale  systems. 

A  third  line  of  work  has  been  followed  by  Campbell  [11]  who  solved 
(1.2)  using  the  Drazin  inverse.  The  Drazin  inverse  is  a  generalized  matrix 
inverse,  closely  related  to  the  eigenspaces  of  a  matrix.  The  solution  of 
(1.2)  can  be  written  in  terms  of  E,  A,  and  B  and  their  Drazin  inverses.  Such 
an  approach  seems  less  traditional  than  that  of  the  others  insofar  as  it 
completely  avoids  explicit  use  of  the  theory  of  matrix  pencils. 

The  motivation  for  descriptor  variable  theory  for  the  most  part  is 
that,  in  choosing  variables  in  a  physical  system  in  a  natural  way,  one  is  not 
always  guaranteed  that  the  resulting  mathematical  model  will  be  in  state 
variable  form.  In  many  instances  the  most  natural  choice  of  system  variables 
may  be  a  non-minimal  set  leading  to  a  system  model  (1.2)  with  E  singular. 
Luenberger  gives  several  examples  of  such  systems  in  [8]  and  [9]. 

In  other  situations  it  may  not  even  be  possible  to  write  a  state 
equation  under  any  choice  of  system  variables.  For  example.  Figure  1.1  shows 
a  simple  electrical  network  with  one  energy  storage  element.  A  convenient 
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choice  of  variables  might  be  as  labeled.  However,  the  loop  equations  are 

Xl“~Xl~X2  (1.4) 

0  =  x^  +  u 

which  cannot  be  manipulated  into  state  space  form.  It  will  be  shown  in 
Chapter  2  that  even  though  the  circuit  contains  an  energy  storage  element,  the 
system  order  is  zero.  (1.4)  is  an  almost  trivial  example  of  one  sort  of 
problem  that  can  arise  in  system  modeling.  However,  it  is  clear  that  in  more 
complex  situations  when  one's  intuition  may  not  be  readily  applied,  a 
systematic  approach  to  such  problems  is  essential. 

We  now  turn  to  singular  perturbation  theory.  It  is  difficult  to 
state  precisely  what  constitutes  a  singular  perturbation  problem.  Singularly 
perturbed  systems  are  identified  as  such  if  they  exhibit  certain  characteristic 
features.  First,  there  is  a  dependence  on  a  parameter  of  some  sort,  usually 
real  and  in  some  sense  "small."  Secondly,  a  slight  change  in  the  parameter 
results  in  a  change  in  system  order.  This  often  occurs  when  a  small  real 
parameter  multiplies  a  derivative  of  a  system  variable.  Many  times,  when 
the  parameter  is  set  to  give  a  low-order  system,  a  descriptor  variable  system 
results.  This  fact  is  the  basis  for  the  relationship  between  descriptor 
variable  and  singularly  perturbed  systems. 

Existing  work  in  singular  perturbation  theory  is  extensive.  Surveys 
of  the  subject  include  [12] -[14]. 


1.2.  Contributions  of  this  Thesis 

In  this  section  we  shall  explore  some  of  the  limitations  of 
existing  singular  perturbation  theory  and  general  directions  that  will  be 
■■aken  in  subsequent  chapters  to  overcome  them. 


The  majority  of  control  ori,  -ed  results  in  singular  perturbation 


theory  have  been  derived  for  the  linear  time- invariant  system 

*1  =  AUX1  +  A12X2  +  Blu 
tox?  =  A2^x^  +  A£2*2  +  B2U 


(1.5) 


where  w  is  a  small  real  parameter  and  A22  is  nonsingular.  There  are,  however, 
important  cases  in  which  system  models  yield  singular  k^-  A  simple  example 
is  given  in  Figure  1.2.  The  system  equations  are 


-01X2  ■  +  u. 


(1.6) 


Here  =  Q  and  c^e  resu^Cs  [  15 ]  —  [  1 7 ]  are  not  applicable.  There  are  only 

a  few  results  available  concerning  singular  k^  (see  [20]). 

There  are  many  instances  in  which  systems  may  not  be  conveniently 
modeled  in  the  form  (1.5)  even  if  singularity  of  A 22  is  allowed.  Consider 
the  operational  amplifier  circuit  of  Figure  1.3.  Assuming  an  ideal  amplifier, 
the  system  equations  are 


<JkiXl  +  X2 


wx2  = 


(1.7) 


(1.7)  seems  to  be  the  most  natural  (i.e.  most  intuitively  meaningful)  de¬ 
scription  of  Figure  1.3.  Yet  it  is  not  clear  how  to  transform  (1.7)  into 
the  form  (1.5).  One  might  try  to  diagonalize  the  matrix  of  coefficients  of 
x^  and  but  this  would  lead  to  a  similarity  transformation  which  is 
singular  at  10  =  0.  Thus  system  equivalence  between  (1.7)  and  the  transformed 
system  would  be  lost  at  w=0.  Other  attempts  at  standardizing  the  model 
must  all  lead  to  a  loss  of  the  natural  interpretation  inherent  in  (1.7)  of 
the  parameter  and  system  variables. 
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Actually,  a  more  natural  interpretation  of  the  system  of  Figure  1.3 
would  include  two  parameters,  one  for  each  parasitic  capacitance.  This 
brings  us  to  the  problem  of  generalizing  the  system  parameter.  It  is  easy 
to  think  of  many  examples  where  the  most  intuitively  pleasing  model  formula¬ 
tion  requires  several  parameters.  However,  little  work  has  been  done  in  this 
area.  Most  results  can  be  applied  only  to  systems  of  a  highly  restrictive 
form  [18] , [19 ] . 

Singular  perturbation  theory  is  primarily  a  qualitative  theory. 

Its  purpose  is  to  give  insight  into  the  nature  of  the  perturbations  of  system 
related  quantities  that  occur  as  a  result  of  slight  changes  in  system 
components.  Quantitative  results  are  scarce.  In  fact,  the  existence  of 
practical  numerical  bounds  on  the  variation  of  system  related  quantities 
would  render  most  convergence  results  obsolete  since  bounds  contain  much 
more  information  than  a  simple  statement  of  convergence.  Singular  perturba¬ 
tion  theory  is  mainly  a  means  of  obtaining  insight  into  the  variational 
characteristics  of  a  system  and,  in  particular,  of  obtaining  information 
about  a  high-order  model  by  examining  one  of  lower  order. 

Given  a  particular  physical  system  it  is  clear  that  one  may  mathe¬ 
matically  characterize  it  with  a  large  number  of  parametrically  dependent 
models.  One  task  of  those  who  apply  the  theory  is  that  of  choosing  the  model 
that  supplies  the  most  information.  Since  the  information  to  be  gained  is 
basically  qualitative,  it  makes  the  most  sense  to  choose  system  parameters 
and  variables  which  have  direct  and  intuitively  clear  relationships  with 
physical  quantities  in  the  system.  It  is  true  that  singularly  perturbed 
models  of  an  extraordinary  nature  may  in  some  cases  be  manipulated  into 
something  resembling  a  standard  form.  However,  it  is  unavoidable  that  such 
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manipulations  at  times  must  diminish  the  intuitive  power  of  the  model 
resulting  in  a  loss  of  information  provided  by  the  theory.  It  is  for  this 
reason  that  a  more  general  theory  is  needed. 

The  main  contributions  of  this  work  are  1)  the  reformulation  of 
the  singular  perturbation  problem  to  include  system  models  which  cannot  be 
naturally  put  into  a  standard  form,  2)  the  extension  of  existing  control- 
theoretic  results  of  standard  singular  perturbation  theory  to  the  more  general 
class  of  system  models,  and  3)  the  unification  of  descriptor  variable  and 
singular  perturbation  theory  for  linear,  time-invariant  systems.  It  is  hoped 
that  this  thesis  will  be  viewed  as  a  fresh  look  at  the  singular  perturbation 
problem.  Whenever  possible,  results  are  stated  in  coordinate-free  or  geometric 
terms.  This  is  done  in  order  to  increase  conceptual  clarity.  In  this  way 
the  reader  is  freed  from  the  burden  of  having  to  keep  track  of  changes  of 
coordinates  which  would  appear  in  any  discussion  of  structural  properties  and 
would  necessarily  depend  on  the  perturbation  parameter.  The  value  of  the 
geometric  approach  has  been  established  by  Wonham  [21]  and  others.  Part  of 
the  last  chapter  is  devoted  to  alternative  algebraic  formulations  of  the 
problem.  This  is  intended  to  give  further  insight  into  the  nature  of 
singularly  perturbed  systems. 

It  should  be  stressed  here  that  even  in  the  linear  time-invariant 
case  the  structure  of  generalized  singularly  perturbed  systems  can  be 
extremely  complex.  Many  pertinent  questions  cannot  be  answered  easily. 


1.3.  Chapter  Survey 

The  thesis  is  divided  into  two  main  parts:  ideas  concerning 
descriptor  variable  theory  (Chapters  2  and  3)  and  singular  perturbation  theory 


(Chapters  4-7).  In  Chapter  2  the  decomposition  theory  of  Rosenbrock  [4] 
is  interpreted  in  geometric  terms.  It  is  shown  that  there  is  a  natural 
decomposition  of  the  state  space  into  independent  subspaces  that  span  the 
whole  space,  including  Rosenbrock's  decomposition  on  the  system.  Chapter  3 
carries  the  geometric  approach  further  by  defining  controllability  of  descripto 
systems  and  studying  the  interplay  between  controllability  and  linear  non¬ 
dynamic  feedback.  The  geometric  structure  of  closed-loop  systems  is  studied 
in  the  tradition  of  Wonham. 

In  Chapter  4  we  begin  the  study  of  generalized  singularly  perturbed 
systems  with  the  definition  of  such  systems  and  with  an  extension  of  the 
geometric  decomposition  described  in  Chapter  2  to  a  region  of  the  parameter 
space.  Under  this  decomposition  a  singularly  perturbed  system  consists  of 
two  subsystems  which  are  consistent  with  the  partitioning  of  the  system 
eigenvalues  into  slow  and  fast  modes.  Chapter  5  is  a  study  of  the  variation 
of  the  trajectory  or  solution  of  a  system  under  small  pertutbations .  Some 
conditions  are  given  under  which  a  small  perturbation  in  system  parameters 
results  in  a  small  change  in  the  system  trajectory.  Chapter  6  is  a  study  of 
the  behavior  of  certain  basic  structural  properties  such  as  stability, 
controllability,  and  stabilizability  that  results  from  a  perturbation  of 
the  system. 

The  linear  quadratic  regulator  problem  is  considered  in  Chapter  7. 
Conditions  are  established  under  which  the  optimal  control,  trajectory,  and 
cost  change  only  slightly  as  a  result  of  a  small  perturbation  in  system 
parameters.  Finally,  Chapter  8  contains  some  alternative  ways  of  looking  at 
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descriptor  variable  and  singularly  perturbed  systems.  Previous  results  are 
interpreted  in  new  ways  and  basic  conclusions  of  the  thesis  are  stated. 
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PART  ONE 

DESCRIPTOR  VARIABLE  THEORY 


CHAPTER  2 


GEOMETRIC  DECOMPOSITION 


2.1.  Problem  Formulation 

To  facilitate  the  development  in  later  chapters  and  to  gain  insight 
into  the  structure  of  descriptor  systems  we  shall  be  concerned  in  this  chapter 
with  the  decomposition  of  a  given  system  into  two  subsystems.  The  decomposi¬ 
tion  has  already  been  achieved  in  [4],  for  example,  using  analytic  techniques. 
However,  the  geometric  structure  of  the  decomposition  has  not  been  described 
elsewhere.  We  shall  parallel  the  development  of  the  analytic  decomposition 
for  descriptor  systems  with  geometric  interpretations  given  at  each  step. 

It  will  be  seen  that  the  natural  response  of  one  subsystem  is  a 
linear  combination  of  the  Dirac  delta  and  its  derivatives.  Hence  it  will 
be  called  the  "fast"  subsystem.  The  other  subsystem  will  be  called  "slow" 
since  it  is  a  state  variable  system  with  exponential  natural  response.  It  is 
emphasized  that  the  terms  slow  and  fast  refer  merely  to  the  natural  responses 
of  the  two  subsystems  of  the  original  open  loop  system.  The  character  of 
the  trajectories  of  the  two  subsystems  may  change  drastically  when  feedback 
or  an  external  control  is  applied.  The  terms  "dynamic"  and  "static"  seem 
to  be  preferred  by  some  authors.  However,  since  we  are  ultimately  concerned 
with  singular  perturbation  theory,  the  terms  slow  and  fast  seem  more 
appropriate . 

Let  X  and  U  be  complex  Euclidean  spaces  with  dim  X »  n  and  dimU*m. 

2  2 

X  and  G  have  inner  products  and  norms  related  by  <  x,x>  =  H  xH  and  <  u,u>  *  ||u|| 
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1  2 

for  x£  X,  u€  U .  ’  We  shall  consider  only  linear  maps  associated  with  X  and  U. 

Let  Hom(U,X)  be  the  C-vector  space  of  homomorphisms  from  U  into  X  where  C  is 

the  complex  plane.  We  choose  to  make  the  distinction  between  homomorphisms 

and  matrices  in  keeping  with  the  philosophy  of  coordinate-free  representations. 

Hom(X,X)  is  the  C-algebra  of  endomorphisms  on  X.  Hom(X,X)  has  identity  element 

I.  If  A€  Hom(X,X)  and  S  is  an  A-invariant  subspace  of  X  then  let  A|S€Hom(S,S) 

be  the  restriction  of  A  to  S.  The  identity  map  in  Hom(S,S)  will  also  be 

denoted  by  I.  If  S  is  not  A-invariant  then  a| S €  Hom(S ,X) . 

Let  A^eHomCX.X),  i=0,...,v  and  choose  a  basis  &  of  X.  We  denote 

by  det(AysV +  . . .  +  A^s +  Aq)  the  determinant  of  the  polynomial  matrix 

[a^  ,sV+  •  •  •  +  a}.s  +  a?.  ]  where  Mato  A.  =  [a^.  ]  .  det(A  s V  +  •  •  •  +  A.,  s  +  A  )  is 
1J  1J  ij  »  ic  ij  v  1  o 

independent  of  the  particular  choice  of  )i  . 

We  shall  consider  dynamical  systems  described  by  differential 
equations  of  the  form 

AxV+'*-+Ax  =  BuTr+**-+Bu  (2.1) 

V  O  TT  O 

where  B^eHomCU.X),  i*0,...,ir  and  superscripts  denote  differentiation  with 
respect  to  t.  A^x'1'  and  B^u1  are  to  be  interpreted  in  the  obvious  pointwise 
sense.  The  class  of  admissible  controls  is  taken  to  be  the  set  of  all  general¬ 
ized  functions  (or  distributions;  see  [24]  and  Section  5.1)  with  range  in  U  that 
are  identically  zero  on  (-*>,0).  The  admissible  controls  form  a  C-vector  space 

denoted  by  £  (U) . 

o 

^"Throughout  the  thesis  all  inner  products  are  denoted  by  (  •,*>, 
and  all  norms  by  1  •!!  . 

2 

<•,•>  is  conjugate  symmetric. 
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The  system  equation  (2.1)  can  be  simplified  considerably  by 

defining  Z*XV*  1JT,  *  (X)  and  (Z)  analogous  to  JJ  (U)  ,  and  z&M  (Z) , 

0  0  0  0 

v€^o(U)  according  to 


z  =  (x, . . . ,x 


v-1 


n-1 

.u  ) 


(2.2) 


and 


v  =  u 


(2.3) 

where  x  is  a  solution  of  (2.1)  for  a  given  u.  Then  (2.1)  may  be  rewritten 


Ez  =  Az  +  Bv 


(2.4) 


where  E,A€  Hom(Z ,Z)  and  B€Hom(U,Z).  Applying  a  control  v  to  (2.4)  corresponds 
to  applying  the  irth  indefinite  integral  of  v  to  (2.1).  Thus  any  problem 
involving  (2.1)  can  be  reduced  to  one  in  which  (2.4)  is  considered. 

Henceforth  we  shall  concern  ourselves  with  (2.4)  only.  Of  course, 
the  interesting  case  occurs  when  E  is  singular.  We  shall  always  assume  that 
det(Es-A) *  0. 


2.2.  System  Decomposition 

For  simplicity,  assume  E,AeHom(X,X)  and  BeHora(U,X).  Our  goal 
is  to  decompose  the  system 


Ex  =  Ax  +  Bu  (2.5) 

k 

into  slow  and  fast  parts.  Let  det(Es-A)  =  it  (s-A  )  where  $  ^  0  and  i  t  j 

0  i=l  1  0 

implies  A^  #  A^  .  Define  a(E,A)  =  {A^, . . .  ,A^}  and  let  A6  C-o(E,A) .  Then 

k  n 

det(AE-A)  =i$  tt  (A-A.)  fO  (2.6) 

o  i 


so  AE-A  is  invertible.  Define 
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and 


where 


S  =  ©  Ker((AE-A)_1E-— D  1 


F  =  Ker((AE-A)_1E)n_r 
k 

r  =  Z  n.  =  deg(det(Es-A))  <  n. 
i=l  i 


(2.7) 

(2.8) 

(2.9) 


Clearly,  r<  rank  E  with  r =  n  if  and  only  if  E  is  invertible. 

Theorem  2.1  gives  a  canonical  decomposition  of  X  with  respect  to  the 
pair  (E,A) . 

Theorem  2.1:  1.  S©F=Xwith  dimS=r. 

2.  There  exists  an  invertible  Me  Hom(X,X)  such  that 

a)  S  and  F  are  both  ME-  and  MA-invariant 

b)  me|s*  i,  ma|f=  i 

c)  ME[F  is  nilpotant 

k  n. 

d)  det(Is-MA|  S)  =  tt  (s-A.)  1. 

i=l  1 

Proof  Let  det(Is-( AE-A)  ^E)  =  sn  d  tt  (s-n.)  ^  where  d=  Z  m.  ,  n.^0  for 

i=l  1  i=l  1  1 

i=l , ,6 ,  and  i^j  implies  rn  t  ^  n-d  is  the  multiplicity  of  the  zero  eigen¬ 
value  of  (AE-A)  '''E.  Define 


EL  =  ®  Ker((AE-A)_1E  -  n.I)  1 

1  i=l  1 


m. 


and 


R2  =  Ker((AE-A)_1E)n“d. 


-1. 


Then  ©R2  *  X,  dim  R^=d,  and  R^  and  R2  are  both  (AE-A)  E-invanant 


p.  28. 


JMuch  of  this  proof  was  patterned  after  Gantmacher  [2],  Vol .  1, 


n.  clearly  depends  on  A.  However,  A  is  fixed  so  we  do  not  write 


this  dependence  explicitly  as  n^(A). 
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Let  J1  =  (\E-A)_1E' R1  and  J2=(\E-A)  1EIR2.  Then  det(Is-J  )  * 

=  T  (s-n.)  1  and  J0  is  nilpotent.  Since  (AE-A)  A=  X(AE-A)  iE-I, 
i=l  i  ^ 

(\E-A)  :a| Rx  =  XJ^-I  and  (XE-A)  ^a|R7  =  XJ2-I.  Define  MS  Hom(X,X)  according  to 


Mx 


r  J~lx 


if  x6  R. 


-1 


( X J I)  X  if  x£  R2 


and  XJ2~I  are  invertible  since  has  no  zero  eigenvalues  and  J2  is 
nilpotent.  Let  M=M(\E-A)  Then  R^  and  R2  are  both  ME-  and  MA- invariant 


with 


and 


Also , 


ME|R1  =  M(AE-A)  1E|R1  =  =  I 

MA|R0  =  M(AE-A)_1A|R?  =  (AJ7-I)_1(aJ2-I)  =  I. 


ME|R2  =  (AJ2-I)  LJ2 


which  is  nilpotent  and 


ma!ri  =  J1  =  *r~Ji • 


Next  observe  that 


det(Es-A)  =  d£t^MA) 
det  M 


det (Is-MA | R^) det (ME | R2s-I) 


det  M 


n-d  k  n,- 

But  det(ME|  R  s-I)  =  (-1)  so  det(Is-MA|  R  )  =  it  (s-d.)  ,  Also 

l  1  i=1  l 


■  _i  <5  i  mi 

det  ( Is-MA  R  )  =  (Is  -  ( AI-J.  L)  )  =  it  (s-(A-— ))  \ 
1  1  i=l  ni 


Thus,  if  the  are  indexed  properly,  we  have  5  =  k,  nu=n^,  d=r,  and 

X-—  *  A  so  •  Hence  R^  =  S ,  R2  =  F.  This  completes  the  proof. 

ni  1  i 
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The  construction  of  M  in  the  preceding  theorem  is  important  so  it 
will  be  repeated  here.  Let 


J  =  (.\e-a)'1e|s 

J2  =  (XE-Aj^ElF 

and  let  M6  Hom(X,X)  be  defined  by 


Finally,  let 


Mx  = 


(U2-D 


X 

if 

xg  S 

X 

if 

x6  F 

M  *  M(AE-A)_i 
Lf  =  ME|F,  Ls  *  ma(s. 


Note  that 


rank  =  rank  ME  -  r  =  rank  E-r. 


(2.10) 

(2.11) 


(2.12) 


(2.13) 

(2.14) 


(2.15) 


Formulas  (2.7),  (2.8),  and  (2.10)-(2.14)  constitute  a  family  of 
algorithms  for  decomposing  the  pair  (E,A).  The  family  is  indexed  by  the 
parameter  A  which  ranges  over  C-a(£,A) .  It  is  fortunate  that  all  the 
algorithms  give  the  same  end  result.  The  following  two  lemmas  establish  that 
S,  F,  M,  and  consequently  Lg  and  are  independent  of 

Lemma  2.1:  Let  \.6fl(E,A),  3  be  any  positive  integer,  and  Lg  and  L^  be  gener¬ 
ated  by  using  a  fixed  parameter  A€  C  in  the  algorithm  (2 . 7)-(2 .14) .  Then 

1)  Ker((nE-A)_1E - - —  I)S  =■  Ker(A  .  I-L  )S 

n-A.  i  s 

and  1 

2)  Ker((nE-A)-1E)S  =  Ker  L® 


for  all  n€  C-cr(E,A) . 
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Proof :  1)  First  observe  that 

(nE-A)_1E--A-  I  = - h—  ( nE-A)-1  (  X  .  E-A) 

n-1.  n-x.  1 

so  if  M  is  given  by  (2.7)— (2.13)  (using  X,  not  n)  we  have 
Ker((nE-A)_1E-— ~  I)B  =  Ker ( (nE-A)"1  ( X . E-A) )  6 

T"|— a  .  1 

1 

=  Ker  ( (nE-A)  "Hf  Hu  X.  E-A))  6 

=  Ker((nME-MA)~1(X.ME-MA))B 
1 

=  Ker((nl-L  )_1(X.I-L  ))S 
s  is 

®  Ker((nLf-I)_1(X.Lf-I))8 

=  Ker ( (nl-L  )'1(XJI-L  ))S 
s  is 

since  X.L..-I  is  invertible.  Also 
l  f 

((nl-L  )-1(x  I-L  ))6  =  (nl-L  )~S(X.I-L  )3 
sis  sis 

since  (nl-L  )  ^  and  X . I-L  commute.  Hence 
s  is 

Ker ( (nE-A)~ LE - —  I)6  =  Ker(X.I-L  )3. 

H” i  s 

2)  The  argument  parallels  that  of  part  1). 

Ker((nE-A)-1E)S  -  Ker ( (nME-MA) -1ME) 6 

-  Ker (n I-L  )-3®Ker((nL,-I)-1Lf.)3 

s  ft 

»  Ker  Lf 

since  (nL^-I)  ^  and  commute.  This  completes  the  proof. 

An  obvious  corollary  to  Lemma  2.1  is  that  S  and  F  do  not  depend  on 
X,  The  \-ihdependence  of  M  will  be  shown  next. 
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Lemma  2.2:  Let  J^(\),  J^(\) ,  M(X),  and  M(A)  be  given  by  (2.10)-(2.13)  for 
some  choice  of  X.  For  n€C-a(E,A)  define  J^(n),  ^(n),  M(n) ,  and  M(n)  in 
the  obvious  way.  Then  M(n)  =  M(X)  regardless  of  the  choice  of  n* 

Proof :  We  have 

M<A)EjS  -  I,  M(A)A|F  =  I. 

Define 

Lg(A)  =  M(X)A|S,  Lf (X)  =  M(X)E|F. 

Then 

J^n)  =  (nM(X)E-M(X)A)_IM(X)E|S 

=  (ni-L  (X))"1. 
s 

Similarly, 

J2(n)  -  (nLf (X)  -  l)_1Lf(X) 
so 

M(n)  |S  *  nl-  Lg(\) 

and 

M(n) | F  =  nL  (A)  -  I. 

Hence  . 

M(n)(nM(X)E-M(X)A)_i  -  I 


so 


M(n)  *  M(n)(r|M(X)E- M(  A)  A)_1M(  A) 
-  M(X) 


which  is  the  desired  result. 

This  brings  us  to  the  decomposition  of  the  descriptor  system  (2.5). 
Operating  on  both  sides  of  (2.5)  by  M  yields 

MEx  ”  MAx  +  MBu.  (2,16) 


Define  P  SHom(X,S)  and  Q  £Hom(X,F)  as  the  skew  projection  operators  on  S 
along  t  and  on  F  along  S  respectively.  Let 
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B  =  PM3 
s 

and 

B  =  QMB . 

We  may  rewrite  (2.16)  as 


(2.17) 

(2.18) 


where 


x  =  L  : 

K  +  B 

s  s 

s  : 

L.x,  =  x  _ 

+  B,u 

f  f  f 

f 

x  ■  Px 

and 

(2.19) 

(2.20) 


We  now  have  two  systems  acting  on  independent  state  spaces  with  x=x  +xr. 

s  f 


2.3.  Trajectories  and  Initial  Conditions 

If  T£Hom(S,S)  define  e(T)  :  [0 ,*>)  Hom(S ,S)  according  to 

e(T) (t)  =  etT.  (2.21) 

Then  the  solution  of  (2.19)  is  simply 

x  =  e(L  )Px  +  e(L  )*B  u  (2.22) 

s  s  o  s  s 

where  denotes  convolution.  Defining  solutions  of  (2.20)  is  a  more 
complicated  task. 

Let  q  be  the  index  of  nilpotency  of  L^.  In  [11]  Campbell 
showed  that  for  each  q  times  differentiable  u:  [ 0 ,<*»)  U  there  exists  a 
unique  differentiable  x^  :  [0,,M)'-»-F  satisfying  (2.20).  x^  is  given  by 

xf(t)  *  -qi1L^Bfu1(t).  (2.23) 

r  i»0  f  £ 


Note  that  no  initial  condition  is  specified.  This  is  in  contrast  to  the 
family  of  solutions  of  (2.19),  a  particular  solution  being  singled  out  by 
choosing  an  initial  condition  x^sS. 
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It  has  been  suggested  in  [5]  and  [7]  that  by  allowing  solutions 
of  (2.20)  from  ^(X),  a  distinct  solution  may  be  defined  for  each  choice 
of  initial  condition  x^6F.  The  proposed  solution  is 


.i-lT i 
-  2  6  L,x 


i=l 


ffo 


q-1  i 
-  2  L,B 


i=0 


re 


(2.24) 


where  d''  is  the  jth  derivative  of  the  Dirac  delta.  (2.24)  was  obtained 
by  taking  the  Laplace  transform  of  (2.20).  Although  the  Laplace  transform 
approach  is  quite  formal  and  is  not  very  satisfying  intuitively,  we  shall 
take  (2.24)  to  be  the  solution  of  (2.20).  A  more  intuitively  pleasing 
justification  will  be  given  in  Chapter  5.  We  shall  see  that,  for  any 
singularly  perturbed  system  with  (2.20)  as  its  limiting  descriptor  form,  if 
its  solutions  converge  to  anything  then  they  converge  to  (2.24). 

The  parts  of  (2.22)  and  (2.24)  due  to  initial  conditions  alone 
serve  as  motivation  for  calling  (2.19)  the  slow  subsystem  and  (2.20)  the 
fast  subsystem.  Also,  we  shall  henceforth  call  S  the  slow  subspace  and  F 
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CHAPTER  3 

CONTROLLABILITY  AND  POLE  PLACEMENT 

3.1.  Reachability  in  che  Fast  Subspace 

In  this  chapter  we  begin  by  defining  controllability  for  descriptor 
systems.  Although  the  idea  of  extending  the  usual  state  variable  definition 
of  controllability  to  descriptor  systems  is  a  fairly  obvious  generalization, 
it  has  not  been  proposed  elsewhere. 

Some  authors  have  come  close  to  considering  controllability. 
Rosenbrock  in  his  theory  of  infinite  decoupling  zeros  [4]  considered 
certain  properties  of  descriptor  variable  systems  which  are  related  to 
controllability.  According  to  his  definitions,  a  descriptor  system  has 
infinite  decoupling  zeros  if  and  only  if  the  matrix  [sA-E  ■  B]  loses  rank 
at  s*0.  It  will  be  shown  that  this  condition  is  equivalent  to  uncontrol¬ 
lability  of  the  fast  subsystem  as  we  shall  define  it  in  Definition  3.1. 
Rosenbrock's  theory,  however,  does  not  address  the  problems  of  state  reach¬ 
ability  and  of  finding  controls  that  steer  the  trajectory  to  a  specified  state. 

In  [10]  Luenberger  et  al.  defined  the  concept  of  maintainability 
which  is  also  related  to  controllability  as  we  shall  define  it.  Maintain¬ 
ability  guarantees  that  a  solution  exists  to  a  certain  type  of  tracking  problem 
determined  by  the  parameters  of  a  descriptor  system.  It  is  clear  after 
examination  of  the  definitions  that  maintainability  is  not  equivalent  to  our 
forthcoming  definition  of  controllability.  For  our  purposes  it  will  be  con¬ 
venient  to  define  controllability  for  descriptor  systems  in  a  way  more  closely 
analogous  to  the  standard  definition  of  controllability  for  state  variable 
systems . 


-  k 
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In  generalized  function  theory  it  is  often  impossible  to  talk  about 
the  value  of  a  function  at  a  given  point.1  To  avoid  this  problem  let  Cq(U) 
be  the  C-vector  space  of  q  times  continuously  differentiable  mappings  from 
[0 ,°°)  into  U  (using  the  right-hand  derivative  at  0)  and  consider  only 
controls  from  Cq(U).  Then  x^  may  be  identified  with  a  differentiable  ordinary 
function  on  (0,°°),  namely 


q-1  i  i 

-  s  9  u1, 

i=0  f  f 


t>  0 


(3.1) 


and  it  makes  sense  to  say 


xf  (t) 


q- 1  2.  i 
-  I  LXu1  ( t) 
i=0  f  f 


(3.2) 


for  any  t>  0. 

The  definition  of  a  reachable  vector  for  descriptor  variable 
systems  is  highly  analogous  to  that  of  a  reachable  state  for  state  variable 
systems.  Let  $  :  [0,”)  *Cq(U)  xX-*-X  be  given  by 


tL  t  (t-r)L 

<f(t,u,x  )  =  e  S(Px  )  +  /  e 
°  °  0 


c  q-1  X  i 

SB  u(T)dT-  1  L3  u  (t) 
s  i=*0  f  f 


(3.3) 


where  u1(0)  is  the  ith  right-hand  derivative  at  0.  Then  ®(*,u,x  )  agrees  on 

(0,®)  with  the  solution  of  (2.5)  with  control  u  and  initial  condition  xq. 

Definition  3.1:  A  vector  wS  X  is  said  to  be  reachable  from  x  £X  with  respect 
-  o 

to  the  system  (2.5)  at  time  xS  (0,®)  if  there  exists  a  control  uS  Cq(U)  such 
that  $(t ,u,xq) 3  w. 

Clearly,  when  applied  to  the  slow  subsystem  (2.19),  Definition  3.1 
is  equivalent  to  the  usual  notion  of  controllability  from  state  variable  theory. 


1What  is  5(0)? 
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We  now  consider  controllability  of  the  fast  subsystem  (2.20).  Let 


6i ^  =  ImB^  +  L^ImB^  +  •••  +  L^ 


(3.4) 


Theorem  3.1:  Let  we  F.  The  following  statements  are  equivalent  with  regard 
to  the  fast  subsystem  (2.20): 

1)  w€  fij. 

2)  There  exist  xg  (O,00)  and  xfQg  F  such  that  w  is  reachable  from  xfo 
at  time  x. 

3)  w  is  reachable  from  x£q  at  time  x  for  every  xg  (0,“),  xfogF, 

Proof:  It  should  be  noted  at  the  outset  that  the  initial  vector  x„  has  no 
-  fo 

effect  on  the  solution  of  (2.20)  for  t>0.  It  appears  in  statements  2)  and 
3)  merely  to  preserve  the  form  of  the  analogous  theorem  from  state  variable 
theory . 

First  we  show  that  1)  implies  3).  Let  gF  and  x£  (0,°°)  be 

to 

given.  Since  wgdL.,  w  =  w  +  •  •  •  +  w  ,  with  w.  6  L^ImB,. ,  i  =  0  , .  .  .  ,q-l.  Since 
f  o  q-1  iff  ^ 

L^ImB^  =  Im(L^B^)  ,  there  are  u.gll  satisfying  L^B^u^-w^,  i  =  Q,...,q-l. 
Define  u€  C^(U)  according  to 


u ( t )  *  UQ  +  (t-x)u1  +  u 2  +  • 


•  +  1^)—  U  . 

(q-1)!  Vr 


Then 


«(x,u,x  ) 


q-l  t  i 
-  Z  L*Bfu  (t) 
i=0  f  f 

q-i  ± 
i*0  f  £  1 


=  w. 


Obviously,  3)  implies  2)  so  it  remains  to  show  that  2)  implies  1) . 
This  follows  almost  trivially  since  L^B^u^(x) S  L^ImB^.  Inspection  of  the 
definitions  of  i  and  gives  the  desired  result. 
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3.2.  Controllability  of  Descriptor  Variable  Systems 

We  now  consider  reachability  and  define  controllability  for 

2 

descriptor  systems.  Let  be  the  controllable  subspace  of  the  slow  sub¬ 
system  (2.19) . 

Lemma  3.1:  Let  p  be  a  nonnegative  integer,  ;S  (0,“),  and 


Then 


t  «  tL  tL*  , 

Jr  =  /  t2pe  SB  B*e  Sdt. 
Tp  0 


ImJV  =  ft  • 
xp  s 


Proof:  Let  x€  Ker.ty'  for  some  x,p.  Then 

-  xp 


x  tL*  0  x  tL  tL* 

/  ||  tPB*e  Sx||  2  dt  =  /  <  x,t2pe  SBB*e  Sx>  dt 

o  s  0 

t  .  tL  tL* 

=  <  x,  /  t~Pe  SB  B*e  x  dt  > 
0 

-  l  x,^Tpx> 


=  0 


tL 


so 


tPB*e  sx  =  0  for  all  t€  [0,x].  Right-hand  differentiation  p+r-1  times  at 


c  *  0  gives 


B*L*  x  *  0,  i *  0, . . . ,r-l. 
s  s 


(3.5) 


Hence 


x€  1"i~1Ker (B*L*i)  *  ^ima^B  V 
i*0  s  s  i*0  s  s 

r-1  i  i 

=  (  E  Im(L  B  )) 
i-0  s  s 


(3.6) 


2See  [21]. 

2  *denotes  the  adjoint  operator. 

^This  is  an  adaptatioii  of  a  proof  given  on  pp.  35-36  of  [21]. 


27 


so  Kerjy  CR  and 
tp  s 


R  c  KerJ^1 

s  Tp 


Im>*  =  Imjrh 

tp  ip 


since  is  Hermitian. 

■*'  xp 

To  show  the  reverse  inclusion  let  xeR  .  Then  (3.6)  holds  and 

s 

hence  (3.5)  does  also.  Recall  that  there  exist  y  :  [0,t]-*C,  i  =  0,...,r-l 


such  that 


tL  1 
e  3  =  Y  ( t)  I  +  y.(t)L  +  •••  +  Y  ,(t)Lj 
o  Is  r-i  s 


for  all  t €E  [0,t],  Thus 


tPB*e  Sx  =  ri1tPY. (t)B*L*1x  =  0 
s  •  n  l  s  s 


and  Jfy  x=0  so 


x€  IraJ^1  *  Ker^  . 
xp  xp 


Hence  R^C  Ker^  and  Imjf'  CR  .  This  completes  the  proof, 
s  xp  xp  s 


Let  R  =  R  3?R^.  The  next  result  justifies  calling  R  the  controllable 
s  t  - 


subspace . 


Theorem  3.2:  Let  w€  X.  The  following  statements  are  equivalent  with  regard 
to  the  system  (2.5): 

1)  w6R. 

2)  There  exists  xs  (0,®)  and  x  €R  SF  such  that  w  is  reachable  from 

os 

x  at  time  x. 
o 

3)  w  is  reachable  from  xq  at  time  x  for  every  x6  (0,°°),  xQSjis«F. 
Proof :  To  show  1)  implies  3)  let  x  and  xq  be  given  and  let  w  =  ws+w^, 
wg€Rg,  w^€R^.  Choose  u^,  i*0,...,q-l  so  that 

q-1  i 

-  1  LcB,u.  ■  w, 
i-0  f  f  1  f 


28 


and  define  u _ €  CH(U)  by 


( t—T  )  ( t”*T  )  ^ 

uf  (C)  =  uo  +  (c-t)u1  +  u2  +  •  •  •  +  uq_r 


T  (x-t)L 

♦  -  /  e 
0 


B  u ,(t) dt  S 
s  f  s 


and  choose  anv  position  integer  p>-ir— .  Let  x  =  x  +  x,  ,  x  e  <R  ,  x,  6  F. 

J  v  s  2  o  so  fo  so  s  fo 

T  Li 

Since  (R  is  L  -invariant,  w  -e  Sx  ~4>€iR  and  from  Lemma  3.1  there  exists 
s  s  s  so  s 

z€  S  with 


z  =  w  -  e  x  - 
^  Ip  s  so 


Define 

,n  (T-t)L* 

u(t)  =  (x-t)  pB*e  z  +  u  (t). 

s  f 

Then  u1(x)  =  u^x)*^,  i=0,...,q-l  and 

1  i 

-  1  LcBcu  (x)  =  w,. 
i=0  f  f  f 


x  (x-t)L 
t  .  s 

1  e 
0 


x  2p  (x-t)L  ( x— t) L* 

B  u(t)dt  =  /  (x-t)  e  B  B*e  z  dt  +  $ 

3  0  3  3 

=  JH-Tpz  +  $ 


and  hence 


$(x,u,x  )  ■  w  +  w.  =  w. 
o  s  f 

2)  follows  from  3)  trivially.  Inspection  of  the  definition  of  <5> 
gives  that  2)  implies  1)  and  the  proof  is  complete. 

For  obvious  reasons,  if  iR =  X  the  system  (2.5)  is  said  to  be 
completely  controllable  and  (E,A,B)  is  called  a  controllable  triple. 


Implicit  in  the  notation  (E,A,B)  is  the  assumption  that  det(Es-A)*  0. 
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Recall  chat  an  eigenvalue  of  Lg  is  called  a  controllable  mode  if  the  eigen- 

space  of  X.  is  reachable,  i.e. 

1 

Ker(.\.I-Ls)  1  c  «g.  (3.7) 

In  [3]  and  [5]  tests  were  established  to  check  for  the  existence 
of  input  decoupling  zeros.  Although  they  were  not  originally  intended  to 
pertain  to  our  concept  of  controllability,  these  tests  are  related  to  our 
definitions.  We  shall  interpret  them  geometrically. 

Theorem  3.3:  1)  An  eigenvalue  A  eo(E,A)  is  a  controllable  mode  of  the  slow 

subsystem  (2.19)  if  and  only  if 


Proof :  1)  Let  M  be  given  by  (2.13).  Then 

M(In(X.E-A) +  ImB)  =  Im(A  ME-MA)  +  Im(MB) 

*  Im(A  I-L  )  +  Im(A.L,-I)  +  Im(HB) 

L  S  X  L 

*  (Im(  A  .  I-L  )  +  ImB  )«F 

is  s 

since  A^L^-I  is  invertible.  From  state  variable  theory, 

Im( A . I-L  )  +  ImB  =  S 
is  s 

if  and  only  if  A is  a  controllable  mode.  Since  M  is  invertible  the  result 
follows . 

2)  From  (3.4)  it  follows  that  the  fast  subsystem  is  completely 
controllable  if  and  only  if  (L  ,B  )  is  a  controllable  pair  or  equivalently, 

Im(AI-L^)  +  ImB^  =  F 
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for  all  'cC.  This  certainly  holds  for  '^0.  Hence 

M(ImE  +  ImB)  =  Im(ME)  +  Im(MB) 

=  S  © (ImL. + ImB^) 

gives  the  desired  result. 

3.3.  The  Effects  of  Linear  Feedback 

Mow  that  controllability  has  been  defined,  we  can  investigate  its 
bearing  on  problems  associated  with  descriptor  systems.  As  we  shall  see,  in 
some  situations  it  is  necessary  to  allow  controls  in  ^q(U)  -  (U) .  At  first 

it  may  seem  strange  that  controllability  is  a  useful  concept  in  such  cases 
since  it  was  defined  entirely  in  terms  of  C^(U)  controls.  However,  closer 
inspection  reveals  that  controllability  is  essentially  a  structural  property 
independent  of  the  types  of  control  driving  the  system. 

Suppose  one  were  to  apply  to  the  system  (2.5)  a  feedback  control  law 

u(t)  «*  Kx(t)  +  v(t)  (3.8) 

where  K€Hom(X,U)  ano  v€.£  (U)  .  The  system  would  then  be  of  the  form 

o 

Ex  =  ( A  +  BK) x  +  Bv.  (3.9) 

However,  it  is  easy  to  construct  examples  where  the  condition 

det ( Es-A-BK)  *  0  (3.10) 

is  violated.  Since  we  do  not  know  how  to  deal  with  such  systems  theoretically, 
only  those  K  satisfying  (3.10)  will  be  considered. 

Our  first  task  is  to  establish  relationships  between  the  structures 


of  the  open  and  closed  loop  systems  (2.5)  and  (3.9).  One  easy  result 
concerns  systems  satisfying 
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rank  E  =  r 


(3.11) 


where  r  is  defined  in  (2.9).  (2.15)  shows  that  (3.11)  is  equivalent  to 

L^=0.  If  (3.11)  holds  then 

dim(Ker  E)  =  n-r  =dimF.  (3.12) 


From  (2.8),  the  definition  of  F,  Ker  ECF.  Hence  (3.11)  implies  F  ■  Ker  : 
so  the  fast  subspace  of  any  system  satisfying  (3.11)  is  invariant  to  linear 
feedback. 

A  more  difficult  result  says  that  the  feedback  invariance  of  the 
controllable  subspace  in  state  variable  theory  can  be  extended  to  descriptor 
systems.  To  prove  this  we  shall  need  a  lemma. 

Lemma  3.2:  Let  ( E , A , B )  have  controllable  subspace  <R.  Then  the  pair  ((\E-A)1E, 
(\E-A)_1B)  has  controllable  subspace  (in  the  state  variable  sense)  for  any 
\€  C-c(E,A) . 


Proof :  Choose  X^o(E,A)  and  observe  that 

R  +  (AI-L  )R  C  R  +  L  R 
s  s 

for  any  subspace  R  of  S.  If  xSR+LgR  then  there  exist  v,  zgR  such  that 
x=v+Lsz.  Let  z  = -z  and  y  =  y+Az.  Then  z,yeR  and 

x  =  y  +  (U-L  )z<=  R  +  ( \  I  —  L  )R 
s  s 

so 

R  +  L  R  =  R  +  (AI-L  ) R. 
s  s 


Assume  that 


R  +  L  R  + 
s 


+  L  R  =  R  +  (\I-L  ) R  + 
s  s 


+  (U-L  )  R 
s 


for  some  k.  Then 
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R  +  L  R  + 
s 


+ 


R  +  L  R  +  • • •  +  LkR 
s  s 

+  L  (R+L  R+  •  •  ■  +  LkR) 
s  s  s 

R  +  L  R  +  ■  •  •  +  LkR 
s  s 

+  (XI-L  )  (R+  L  R+  •  .  .  +  LkR) 
s  s  s 

k+1 

(XI-L  )R+  •  •  +  (aI-L  )  R. 
s  s 


Setting  R=ImBs>  it  follows  that 

<R  =  ImB  +  (XI-L  ) ImB  +  •••  +  (XI-L  )n_1ImB  . 
s  s  s  s  s  s 

Since  iR  is  L  -invariant  it  is  also  (AI-L  )  ^-invariant  so 
s  s  s 

<R  =  (XI-L  )_1ImB  +  •••  +  (XI-L  )“nimB  . 

S  3  g  g  S 

We  shall  now  prove  that 

flf  =  (XLf-I)~1ImBf+  ((ALf-I)'1Lf)(XLf-I)~1IinSf  +  ••• + 

+  ••  .  +  ((XLf-I)"ILf)I1_1(XLf-I)'1ImBf. 

For  X  =  0  this  is  obvious  so  assume  that  X  £  0  and  let  R  be  any  subspace  of  F. 
Clearly , 

Lj_1R+  (XLf-I)Lj~2R+  (XLf-I)2L""3R+  •••+  (XLf-I)n-1RC  R  +  LfR  +  ■  •  •  +  l"_1R. 


Let  0<k<n-l  and  consider  the  nxn  matrix  T  = 

t 


t 


ij 


(  1  )(-l)J  X1_j 
0 


[  t .  .  ]  where 
if  i>  j 
if  i<  j 


V 


T  is  invertible  so  there  exist  a  ,  i =  0, 


n-1 

L 

i=j 


t .  .  'j. 
ij 


i 


0  if 

1  if 


. . ,n-l  satisfying 
j  #  n-k-1 

j  =  n-k-1 


mi  _  ,iw  n+j-i-1 

2;  a  (XL  -1)  L  =  I  l  a  (  )  (-1)  JAJL 

l-O  1  f  f  i=0  j=*0  1  J  f 

HZ1  .  i  .  ,  .  j  A  i-j  n-j-l 

j  ssQ  i=j  i  i“j  f 


If  xG  R  then 


k.  H-l  n-i-1 

L*x  =  I  (\LC-I)  Lc  (a.x) 


eL“_iR  +  (\Lf-I)L£'2R  +  •••  +  ( XL f  —  I ) n  LR 

and 

R  +  LfR+  •••  +  L£_1R  =>  l”_1R  +  (AL£-I)L£“2R  +  •••  +  ( ALf-I)n_1R. 

Setting  R=  ImB£  and  observing  that  (R£  is  (A  L£-I)  "''-invariant  gives 
Slf  *  (ALf-I)~1ImBf  +  (ALf-I)~2LfIm,&f  +  •  •  *  +  (XLf-I>'nL£_1ImBf . 

Note  that 

Lf(ALf-I)  =  ( ALf-I)Lf 

so  left  and  right  multiplication  by  (AL£-I)  2  shows  that  (AL£-I)  2  and  L£ 
commute.  This  establishes  the  desired  expression  for  ft £. 

Finally,  note  that 

(AE-A)-1E|S  *  ( AI-L  )_1 
S 

( AE-A) _1E  j  F  =  (ALf-I)_1Lf 

P ( AE-A) _1B  =  (AI-L  )"1B 

s  s 

QOE-A)“1B  =  ( ALf-I) _13f 

so  the  pairs  ((AE-A)  2E|S,  P(\E-A)  2B)  and  ((\E-h)  2E|  F,  Q(AE-A)  ^B)  have 


controllable  subspaces  iR  and  (R.  respectively.  Since  (AE-A)  E  |  S  and 
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( \ E— A)  1E|  F  have  disjoint  spectra,  ((XE-A)  ^E,  (XE-A)  1B)  has  controllable 
subspace  Rg  3R^  =R  .  This  completes  the  proof. 

Theorem  3.4:  The  triple  (E,A+BK,B)  has  controllable  subspace  R  for  all 
K€  Hom(X,U)  that  satisfy  (3.10). 

Proof:  Let  R„  be  the  controllable  subspace  of  (E,A+BK,B).  For  any  subspace 

-  rv 

R  of  X  and  any  X €  C- (a (E , A) n a (E , A+BK) )  we  have 

(XE-A-BK) (XE-A)-1 (ImB  +  R)  =  (I-BK(XE-A)-1) (ImB +  R) C  ImB +  R 

and  -1  -1 

(XE-A)  (ImB  +  R)  =  (XE-A-BK)  (ImB  +  R) 

since  XE-A-BK  is  invertible.  Applying  Lemma  3.2  gives 

R  =  (XE-A)-1 (ImB  +  E( XE-A)-1 (ImB  +  E(XE-A)-1( • • • (ImB) •••))) 

=  ( XE-A-BK) -1( ImB +  E (XE-A-BK)-1 (ImB +  E (XE-A-BK) -1 ( • • • (ImB) •••))) 

=  Rjr 

and  the  proof  is  complete. 

Henceforth,  for  notational  simplicity,  we  shall  denote  the  relevant 
subspaces  and  operators  of  the  closed  loop  system  (3.9)  by  SR,  FR,  Mj^,  Lgj., 

Lf K*  etc' 


3.4.  Slow  Feedback 

Besides  feedback  invariance  of  the  controller  subspace,  there 
does  not  appear  to  be  much  that  can  be  said  in  general  r  '  iting  structural 
properties  of  open  loop  descriptor  systems  to  those  of  closed  loop  ones. 
Fortunately,  the  pole  placement  problem  can  be  dealt  with  by  feeding  back 
the  slow  and  fast  trajectories  separately.  The  induced  structural  changes 


can  then  be  more  easily  characterized.  In  this  section  we  consider  feedback 
in  the  slow  subsystem.  That  is,  we  apply  a  control  u=Kx+v  with 


Ker  K  T>  F  . 


Let 


Ks  =  K|  S 


and  let  a  =  (e,  , . . .  ,e  )  and  )i  ,  =  (e  , ,  , . 

s  i  r  f  r+1 


.,e  )  be  bases  of  S  and  F 
n 


respectively.  If  *  (e,  , . . .  ,e  )  then 

1  n 


Mat  (MEs-MA-MBK)  = 

si 


MatR  (Is-L  -B  K  ) 
»  s  s  s 


Matrt  *  (-BfKs) 
f  s 


Mat^  (Lfs-I) 


(3.13) 


(3.14) 


(3.15) 


Clearly,  the  eigenvalues  of  the  closed  loop  system  (3.9)  are  those  of  the 

operator  L  +  B  K  .  Hence  we  have  the  following  extension  of  a  well  known 
s  s  s 

result  from  state  variable  theory. 

Theorem  3.5:  An  eigenvalue  X^  of  the  descriptor  variable  system  (2.5)  can 

be  shifted  arbitrarily  by  applying  slow  feedback  if  and  only  if  the  eigenspace 

of  X.  is  contained  in  ft  . 

1  s 

From  (3.15)  it  follows  that  the  dimensions  of  the  slow  and  fast  sub¬ 
spaces  do  not  change  when  feedback  is  applied.  In  fact,  the  next  result  says 
that  the  fast  subspace  and  fast  subsystem  are  essentially  unchanged  by  slow 
feedback. 

Theorem  3.6:  If  K  satisfies  (3.10)  and  (3.13)  then 

FK=  f 
LfK  “  Lf 

and  flfK  =  ftf . 


Proof :  Choose  X  €  C-(a(E,A) ^  a(E,A+  BK) )  and  observe  that 
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Mat  ((,\E-A-BK)_1E) 


Mat  .  (\I-L  -B  K  ) 

ii  s  s  s  s 

Mat  ( \L_-I)_1B^K  OI-L  -B  K  ) 
ii  Ji  f  f  s  s  s  s 

t  s 


0 


Mat,- 


(3Lf-I) 


Clearly , 


F  =  Ker((AE-A-BK)“1E)n_r  =  F. 
K 


Also,  (3.13)  implies 

(\E-A-BK)“1e1f  =  (AE-A)-1E|F 

—  —1  —1  — 

mJf  =  (-U2K-I)  =  (XJ2-I)  =  M|F 

and  _ 

LfK  =  EiF  =  M(AE-A)  E j  F  =  Lf. 


Finally,  ^fK  =  ^^nEK,  FK=F,  and  Theorem  3.4  together  imply  (R^K  =  (Rf  so  the 
proof  is  complete. 

It  is  easy  to  construct  examples  where  B^^B^..  Nevertheless,  as 
we  have  just  seen,  the  open  and  closed  loop  systems  have  the  same  fast 
controllable  subspace. 


3.5.  Fast  Feedback 

Consider  the  control  law  u=Kx+v  with 

Ker  K  3  S 


and  let 

Kf  *  K|F. 


(3.16) 

(3.17) 


Then 
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1 

I 

I 

I 


Mat  (MEs-MA-MBK) 


Mat^  (Is-Lg) 
s 

0 


Mat 

Mat , 
a 


ii  ii  ,  St 
s  f 

<LfS-I-BfKt) 


(3.18) 


Clearly,  the  eigenvalues  of  the  open  loop  system  are  also  eigenvalues  of  the 
closed  loop  system.  But  det (L^s-I-B^K^)  in  general  is  not  a  constant  poly¬ 
nomial  so  fast  feedback,  may  induce  additional  modes  in  the  system. 

If  some  of  the  roots  of  det(L^s-I-BfK^)  are  also  eigenvalues  of  Lg 
then  it  is  difficult  to  find  a  relationship  between  the  open  loop  and  closed 
loop  eigenspace  structures.  However,  this  can  be  easily  voided  as  we  shall 
now  see. 

Let  X  denote  the  subspace  of  Hom(X,U)  consisting  of  all  T 
satisfying  KerTCS  and  let  \  consist  of  all  T  satisfying  det  (L^s-I-B  )  *  0 
and  such  that 

c  (Lf,I  +  BfTf)  na(E,A)  =  >  (3.19) 

where  =  l| F .  Let 

II  Til  *  sup{IITxll  |xe  X, Hxll  =  1}.  (3.20) 

Proposition  3.1:  T  is  open  relative  to  X  . 

Proof :  The  proof  will  be  postponed  until  Chapter  4.  See  the  discussion 
following  Lemma  4.3. 

Corollary :  There  exists  e>0  such  that  TS?C  and  IITII  <e  together  imply  that 
T  €  I  . 

Prop f :  Obviously  Osl  so  the  result  follows  immediately  from  the  proposition. 

We  next  establish  a  threefold  decomposition  of  the  closed  loop 
system.  Let 


*v-. 
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h  8< 

det(Lcs-I-BrK,)  =  tt  (s-S.) 
f  ft  °i=l  1 


(3.21) 


where  ?  VO  and  i  f  1  implies  3.^3..  Henceforth  we  assume  that  KeT  so  3.^  A. 

o  i  j  1  j 

for  all  i,j.  Let  AG  C-o(E,A+BK)  and  define 


h  _1  1  §i 

D„  =  0  Ker ( ( XE-A-BK)  XE  -  I)  \ 


i=l 


X-6. 

i 


(3.22) 


Lemma  2.1  guarantees  that  is  independent  of  X. 

1)  sk=s©dk. 

2)  S  and  are  both  M^E-  and  M^A+BK) -invariant  with 


Theorem  3.8:  1)  S„  =S©D„. 

r  K.  K. 


M^A+BK)  |  S  -  L  . 

3)  «SK-«,*V 

Proof :  1)  Let  AG  C-a(E,A+BK)  and  define 

_  t  _i  i  n-t 

S  *  8  Ker ( (AE-A-BK)  E  -  -—-)  . 

i=l  A~ 

Then  S^=S^D^.  Let  x€  X  with  x=!x^+X2,  x^GS,  and  Then 

( XE-A-BK) _1Ex  =  ( XI-L  )-1x1  +  (XI-L  )~1B  Kr (XL,-I-BcKc)_1Lcx- 

s  1  s  str  rr  r  z 

+  (XLf-I-BfKf)_1Lfx2 


and  there  exists  Ng  Hom(F,S)  such  that 

n . 


k  .  .  n .  k  .  .  n . 

tt  (  ( XE-A-BK) _iE  -  -H=r—  I)  1x  =  Nx„  +  it  ((XLC-I-BCKC)“  Lc  -  I)  Xx,  . 

i»l  A~ A..  ^ 


i=l 


f  “f~f'  “f  X-X. 

l 


Since  K€T  ,  -  is  not  an  eigenvalue  of  (XL^-I-B^K^)  ^L^  ^  so 


(“•[■‘•W  le‘  — 


I  is  invertible.  Thus,  if  x£  S  we  have  x? *  0  and 


xGS.  Conversely,  suppose  xe  S.  Then 


See  the  proof  of  Theorem  2.1  for  a  discussion  of  the  eigenvalues 
of  (AE-A)“^E  for  any  E  and  A  with  XE-A  invertible. 


? 

i 

a 


Ic  . 

t  ((XE-A-BK)  E--~-  I)  LX  =  0 
i=l  X~A,- 


and  x€  S  so  S  =  S . 

2)  We  have  immediately  that  S  and  Dj,  are  (XE-A-BK)  ^-invariant, 
J  -  and  M^-invariant ,  and  hence  M^E-invariant .  From 

( XE-A-BK) -1(A+BK)  =  A (XE-A-BK) _1E- I 

M^(A+BK) -invariance  of  S  and  follows.  From  (3.16), 

(XE-A-BK)-1E| S  =  (XE-A)_1E|S 

\ls  ■  jIrI  s  "  JI*  ‘  "Is 


M,, (A  +  BK)  |  S  =  SL  ( XE-A-BK)  A|S 

=  m(ae-a)-1a|s 


=  L  . 
s 


3)  Clearly, 


Im(  B  .  ME-MA-MBK)  +  ItnMB  =  S  (Im(B.Lc-I-BcK,)  +  ImB,)  . 
i  l  r  r  r  f 

For  x€  F  let  x^  =  (gX^-I)  ^x  and  x^^K^x^.  Then 
(BiLf-I-BfKf)x1  +  Bfx2  =  x 


Im(8iLf-I-BfKf)  +  ImBf  =  F. 

Hence,  from  Theorem  3.3,  part  1),  D,, C  ift .  Also  D„CS,.  so 

N  A  l\ 

DK  c  s  n  SK  ■  flsK- 
Furthermore,  by  Theorem  3.4, 


<R  »«nsc«  ns 

s  K  K  sK 
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so 


flsK  D  «s®V 


To  prove  the  converse  let 


x€flsfC  =  «Kn  sk  =  a  n  (S  ®Dr)  . 

Then  xSfl  and  there  exist  y€S,  z£D  such  that  x=y+z.  But  y  =  x-z6fl  +  D  =<R 

K.  K 

SO 


and 


xS  (fins)  0D„  =  <R  0D„ 

K.  S  K 


«sKC«s®dk. 


This  completes  the  proof. 

Theorem  3.8  is  analogous  to  Theorem  3.6.  It  states  that  the  closed 

loop  system  consists  of  three  subsystems:  one  acting  on  the  open  loop  slow 

subspace  S  with  eigenvalues  X^  and  controllable  subspace  (Rs ,  one  which  is 

completely  controllable  acting  on  D  with  the  induced  eigenvalues  and  a 

K  1 

fast  subsystem  acting  on  F  .  Although  there  is  considerable  structural 

K. 

reshuffling,  the  controllable  subspace  of  the  overall  system  <R  remains 
unchanged. 


If  rank  E=r  then  we  can  go  even  further.  We  already  know  that  in 
this  case  the  fast  subspace  does  not  change  when  feedback  is  applied.  Since 
SC  S„  and  S  =  X  it  follows  that  S  *  S  .  The  next  result  says  that  not  only 
are  the  slow  and  fast  subspaces  unchanged,  but  the  entire  system  is  essentially 
unaffected.  This  would  seem  to  indicate  that  applying  fast  feedback  to  such 
a  system  is  pointless.  Recall  that  rank  E=r  implies  ■  0  so  I  +  B^K^  must 
be  invertible  for  det (Lcs-I-B^K^)  *  D  to  hold. 

Theorem  3.9:  If  rank  E=r  and  K  satisfies  (3.10)  and  (3.16)  then  L  „=L_, 

SK.  S 

hit*0’  BsK-Bs(I-Kf<1+BfK£)'1V'  and 
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Proof  :  Since  F=KerE,  =  >I^E |  F  =  0 .  Also,  if  \S  C-a(E,A+BK)  then 
J.v  =  (XE-A-BK)_1e|S  =  (AE-A)~1E | S  =  (AI-L  )_1 

IK  S 

and  , 

LsK  =  \(XE~A-BK)~  ^A+ BK)  Is 

=  j~*(ae-a)-1a|s 

=  L  • 
s 

Next  we  have 

™K  ■  4p 

so  -1  -1 

PM  B  =  J^PUME-MA-MBK)  MB 
K  J.K 

=  Bg(I-Kf(I  +  BfKf)_1Bf). 

Finally , 

nc  *  <aj2  k-i)_1q 

so  -l  -l 

QMrB  =  (XJ^-I)  ■‘■Q(AME-MA-MBK)  mb 

-  -(XJ2K-I)-1(I  +  BfKf)_LBf. 

But  , 

J2K  =  (  E-A-BK)"  E ! F  =  0 


which  gives  the  desired  result. 

To  conclude  this  chapter  we  consider  the  problem  of  eliminating  the 
impulsive  portion  of  the  fast  trajectory  (2.24)  by  applying  fast  feedback. 

The  result  that  we  shall  obtain  says  that  it  is  possible  to  eliminate 
impulsive  behavior  if  and  only  if  the  fast  subsystem  (except  for  the  part 
acting  on  Ker  L^)  is  controllable.  First  we  need  a  lemma. 

Lemma  3.3:  Let  Y  and  Z  be  finite-dimensional  C-vector  spaces  with  dimY=dimZ. 
Let  NSHom(Y,Z)  and  GeHom(U,Z).  There  exists  H€  Hom(Y,U)  such  that  N  +  GH  is 
invertible  if  and  only  if  ImN+ImG=Z. 


Proof:  If  N  is  invertible  then  the  result  is  obvious.  Let  N  be  noninvertible 


and  choose  bases  of  Y,  Z,  and  U.  The  existence  of  an  appropriate  H  is 
equivalent  to  controllability  of  the  mode  0  of  the  pair  (MatN,  MatG)  which 
is  equivalent  to 

ImN  +  ImG  =  Z 


so  the  proof  is  complete. 

Let  F/KerL^  be  the  quotient  space  of  F  modulo  KerL^, 

Wg  Hom(F,F/Ker  L the  canonical  surjection,  and  L^  the  induced  map  of  L^. 
Proposition  3.2:  Lf  is  nilpotent  with  index  of  nilpotency  q-1. 

Proof:  L^  is  uniquely  defined  by  WL^*L^W.  Assume 


Then 


WLp  =  L»W. 
f  f 


WL^+1  *  L^WLf 


(3.23) 


lP+1w. 


Hence  (3.23)  holds  for  p=  1,2,3,... 

Next,  note  that  if  x€  then  there  exists  y£  F  with  x=L^  ^y. 

Thus  L^x  =  L^y  *  0  and 

ImL^-1  C  Ker  Lf  =  Ker  W 


so 


L^_1W  =  WL^'1  =  0 


and,  since  W  is  a  surjection,  L 


q-i 


0. 


r  q-1 


Thus 


On  the  other  hand,  there  exists  ye  F  with  Lj;  y  f4  0  so 
ImLg-2  £  Ker  Lf. 


L^'2W  =■  WL^"2  i  0 


and  L^_2?t0.  This  completes  the  proof. 
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Letting  B^  =  WB^  and  x^(t)  = Wx^(t) ,  we  may  define  the  quotient  system 

+  B^u.  (3.24) 

Since  is  nilpotent,  (3.24)  has  similar  structural  properties  to  those  of 
the  fast  subsystem  (2.20).  In  particular,  from  Theorem  3.3,  part  2)  it 
follows  that  (3.24)  is  completely  controllable  if  and  only  if 
ImL^ +  ImBf =  F/Ker  L ^ .  We  shall  make  use  of  this  fact  in  the  final  theorem 
of  this  chapter.  Note  that  no  5-functions  are  present  in  (2.24)  if  and  only 
if  *  0. 

Theorem  3.10:  The  following  statements  are  equivalent: 

1)  There  exists  KS  Hom(X,U)  satisfying  (3.10)  such  that  L  =0. 

tK. 

2)  There  exists  K€  Hom(X,U)  satisfying  (3.10)  and  (3.16)  such  that 
LfK=°- 

3)  ImLc  +  ImB.  +  KerL, *  F. 

ft  t 

4)  The  quotient  system  (3.24)  is  completely  controllable. 

Proof :  Choosing  bases  of  S  and  F,  it  is  clear  from  the  matrix  representation 
of  Ex-MA-MBK  that,  for  any  K, 

deg(det (Es-A-BK) )  =  r  +  deg(det (L^s-I-B^K^) ) . 

From  (2.15),  L^=0  if  and  only  if 

deg(det(Es-A-BK) )  =  rank  E 


or  equivalently, 


since 


deg(det (LjS-I-B ^K^) )  =  rank 


rank  E  *  r  +  rank  L  ^ . 


Whether  or  not  L  *  0  is  therefore  determined  solely  by  the  action  of  K  on  F. 
t  K. 
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The  behavior  of  K  on  S  is  irrelevant  and  the  equivalence  of  1)  and  2) 
follows . 

Choose  a  basis  =(e  ,...,e  ,e  , . . . ,e  , . . .  ,e  ,...,e  ) 

1  P1  Pl+1  P2  Pd-1+1  Pd 

of  F  so  that  Mato  L^  is  in  Jordan  form  with  d  blocks  of  sizes  p.,,-p.,  Let 
»  f  l+l  l 

Mat  (I  +  B,Kf)  =  [h,.].  A  straightforward  but  notationally  messy  calculation 

Et  1J 

yields  that  the  (rankL^)th  coefficient  of  det(L^s-I-B^K^)  is  just  detQ 
where,  setting  po  =  0,  3=  [9  ]  with 

0 .  .  =  h  , .  . 

ij  pi’pj-l+1 

1)  is  equivalent  to  det  3^0  for  some  K  .  Note  that 


and 


Let 


ImLf  =  span{e_.  |  j  =  pi_1+l,.  .  .  ,pi~l;  i=l,...,d} 


Ker  L£  =  span{ere  r...,e 

1  a- 1 

T  =  span{e  ,e  , . . .  ,e  } . 

P1  P2  Pd 


}. 


Then  dim(  Ker  L^)  =  dim  T  and  1)  is  equivalent  to  the  statement  that 

P__  T  (I  +  B..K..) J Ker  L,  is  invertible  for  some  Kc  where  P„T  .  is  the  skew 
llmLj  t  t  t  t  iimL^ 

projection  operator  on  T  along  IraL^.  Let  V=PTImLf-  Th6n 

V(I  +  BfKf)  |Ker  Lf  -  v|KerLf+  (VBf )  (Kf  |  Ker  L  ) 

and  from  Lemma  3.2  an  appropriate  K^  may  be  found  if  and  only  if 

V(Ker  Lr  +  ImB,)  =  ImV|KerLr  +  ImVB,  =  T. 
ft  f  t 

Hence  we  have  arrived  at  the  equivalence  of  3) . 

Complete  controllability  of  (3.24)  is  equivalent  to 


ImL^  +  IraB^  =*  F/Ker  L^ 


so  the  equivalence  of  3)  and  4)  follows  from  elementary  arguments. 


45 


Since  1)  implies  2)  in  Theorem  3.10,  we  are  guaranteed  that  if  the 
impulsive  behavior  of  the  fast  subsystem  can  be  eliminated  by  feedback  then 
it  can  be  eliminated  by  fast  feedback.  Theorme  3.10  may  be  interpreted  as 
a  pole  concerned  with  shifting  poles  at  infinity  into  the 

finite  portion  of  the  complex  plane.  Theorem  3.8  says  that  the  shifted  poles 
correspond  to  controllable  modes  and  can  thus  be  placed  arbitrarily. 


3.6.  A  Two-Stage  Pole  Placement  Procedure 

The  following  design  procedure  can  be  used  for  pole  placement  in  the 
overall  descriptor  system.  First,  calculate  the  decomposition  for  the  given 
open  lcop  system.  If  the  fast  subsystem  modulo  Ker  is  completely 
controllable  then  any  impuslive  behavior  can  be  eliminated  by  applying  fast 
feedback  as  outlined  in  the  corollary  to  Proposition  3.1  and  Lemma  3.3. 

Second,  calculate  the  decomposition  of  the  closed  loop  system  after 
fast  feedback  has  been  applied  and  shift  the  poles  of  the  slow  subsystem  as 
desired.  The  properties  listed  in  Theorem  3.8  make  calculating  the  decompo¬ 
sition  easier. 


\ 
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PART  TWO 

SINGULAR  PERTURBATION  THEORY 
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CHAPTER  4 

DECOMPOSITION  OF  GENERALIZED  SINGULARLY  PERTURBED  SYSTEMS 

4.1.  preliminaries 

As  shown  in  Chapter  1,  there  is  a  need  for  a  singular  perturbation 
theory  which  not  only  unifies  existing  theories,  but  also  extends  them  to 
a  larger  class  of  systems.  In  this  Chapter,  after  defining  the  generalized 
singularly  perturbed  system,  we  shall  develop  a  geometric  decomposition  of 
the  system  into  slow  and  fast  subsystems.  Such  a  decomposition  will  have 
use  in  Chapter  5  in  the  study  of  the  behavior  of  the  solutions  of  (4.8). 

It  will  also  be  useful  in  studying  the  behavior  of  the  solution  to  the  LQ 
regulator  problem  in  Chapter  7.  For  the  standard  form  (1.5),  approximate 
decompositions  already  exist  (see,  for  example,  [15]).  We  shall  take  a 
somewhat  different  approach  from  what  has  been  done  in  the  past,  extending  the 
geometric  decomposition  developed  in  Chapter  2  to  systems  defined  on  a  para¬ 
meter  space.  The  decomposition  will  be  exact  in  contrast  to  the  approximate 
decoupling  result  of  [15].  In  order  to  develop  the  theory  we  shall  need  certain 
mathematical  concepts.  Consider  the  set  /7<  of  all  subspaces  of  X.  For  any 

R  £/7i  let  PR  £Hom(X,X)  be  the  orthogonal  projection  operator  on  R.  Define 
2 

p:  7T(  -  [0,»)  by 

p(R,T)  =  ||PR  -  PTil.  (4.1) 

It  is  shown  in  [33],  pp.  69-71  that  p  is  a  metric  with  values  in  [0,1]  and 
that  p(R,T)  *  1  if  and  only  if  either  R- fi T  /  0  or  R H Tx  # 0.  From  a 
dimensionality  argument  it  follows  that  p(R,T)  =  1  when  dim  R/dim  T.  p  can 
be  thought  of  as  a  generalization  of  the  angle  between  subspaces. 


I 
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Let  be  a  topological  space  (see,  for  example,  [28])  and  choose 

x>0  (P..  If  R  :  Q - 7/(  is  continuous  at  xq  with  respect  to  p,  it  will  be  useful 

to  construct  a  convergent  basis  for  R„  To  do  this  we  need  a  pair  of  lemmas. 

Lemma  4.1.  Let  y,  e^:Q— X,  i  *  l,.,.,n  be  continuous  at  with  (e^(ou) , » . . , 

en(o>))  a  basis  of  X  for  all  xtfi.  If  o^:  Q-C,  i  =  l,...,n  are  defined  by 

y(x)  =or,  (x)e,  (x)  +  . . .  +  or  (x)e  (x)  then  each  a.  is  continuous  at  x  . 

-'v/l  l  n  n  i  o 

Proof ;  Define  L  :  Q  — Hom(X,X)  according  to 

L(x)ei(xQ)  =  ei(x)  ,  i  =  l,...,n. 

Then  L  is  continuous  at  xq  and  L(x)  is  invertible  for  each  xcQ.  Taking 

the  inverse  of  L(x)  corresponds  to  a  continuous  function  on  the  topological 

subspace  of  invertible  endomorphisms  in  Hom(X,X).  Hence  x~*L(x)  ^  is 

continuous  at  x  .  It  follows  from 
o 

||l(x)  ”1y  (x)  -  L(xo)'ly(xo)P  *  !iL  (uu)  "  L-L  (c«o  )  “  L  H  Sly  (a.)  il 
+  !lL(U)o)"1!|!ly(x)  -  y(uu0)ll 

so  x~*L(x)  y(u»)  is  continuous  at  xq.  Define  :  Q-C  according  to 


L(x)  Ly(x)  *  3.  (x)e.  (x  )  +  . . .  +  P  (x)e  (x ) . 

llo  n  n  o 


Then  each  is  continuous  at  x  .  It  follows  that 
i  o 

y(x)  =>L(x)(@L(uu)e1(xo)  +  ...  +  @n((u)en(uJo)) 

=  0.  (x)e1  (x)  +  . .  •  +  @  (x)e  (x) 
ll  n  n 

so  i*l,««*»n  and  the  proof  is  complete. 

Lemma  4,2.  Let  L  :  Q— Hom(X,X)  be  continuous  at  xq  with  (v^,...,v  )  a  basis 
of  Ker  L(xo).  There  exists  a  neighborhood  V  of  xq  and  maps  e^  :  V  — X, 
i  ■  1, . . . ,p  such  that 
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1)  e.  is  continuous  at  n  ,  i  =  L,...,p. 

L  0 

2)  ei(^o)  =  v.,  i  =  l,...,p. 

3)  (e^(uu) , . . .  ,ep(uj))  is  linearly  independent  VuutV. 

4)  span  {e^(uj)  , . . .  ,ep(uj)  j  3Ker  L(uj)  V  uj  £  V. 

5)  (jj  — span  {e.  (<u) , . . .  ,e  (uo)}  is  continuous  at  ou  . 

ip  o 

Proof:  Choose  a  basis  (b,,....b  )  of  ImL(u)  )  and  define  c.  :  U-X. 

-  1  n-p  o  l  * 

i  =  l,...,n-p  according  to  c^(uu)  =L(ou)  b^.  Then  c^((iu)  Ker  L(ou) x.  Further- 

&  ( 

more,  (b.  , . . .  ,b  )  is  a  basis  of  Ker  L(uu  )  -  so  (c.  (a»  ) , , . .  ,c  (uu  ))  is  a 
’  1  n-p  o  1  o’  *  n-p  o 

*  n-p 

basis  of  ImL(u>o)  =  Ker  L (<JUQ ) .  Since  the  points  in  X  which  correspond 

to  linearly  independent  sets  of  vectors  form  an  open  set,  there  is  a 

neighborhood  of  ibq  throughout  which  (c^(uj)  , . . .  >cn_p0"))  is  linearly 

independent.  Taking  the  adjoint  corresponds  to  a  continuous  function  on 

*  *  *  "  , , 

Hom(X,X)  so  ou-L(ou)  is  continuous  at  ou  .  ilc.(w)  -c.(ou  )i!  £  iiL(w)  -  L(uj  )  , i  j) b . , 

o  i.  i.  o  o  i 

gives  that  is  continuous  at  i=l,...,n-p. 

Applying  the  projection  theorem  (see  [27],  p.  56)  yields  the 

orthogonal  projection  p  .  (ui)c.  (ou)  +  . . .  +  p.  _  (ou)c  _  (tu)  of  v.  on  span 

x  jii"p  n-p  l 

Lc^(su)  ,...  ,cn_p(ou)]  CKer  L(ou)x  where  :V^-C  is  given  by 


<cl(ou),ci(uu)>  ...  <cL(m),c  (u))>  PtL(ou) 


<Cn-p(tu)’Cl(u,)>  *'•  <cn-p(u,)»Cu-p(u,)>  @i,n-p(uj) 


<Vh<u*)> 


<v,,c  (ju)> 

l  n-p 


Clearly,  each  p.  .  is  continuous  at  uu  . 

J.J  o 


Define  e^  on  by 


e  (<u)  -v  -  (0  .  (ou)c  (<ji)  +  ...  +p  (a)c  (uu))  . 

j-  i  xi  l  i,n-p  n-p 


Then  each  e.  is  continuous  at  ou  .  Also, 
i  o  ’ 
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v.  span  lc ,  (xi  )  . .  .  c  (jj  ); 
x  -  1  o  n-p  o  J 

so  e^O^)  =  v^,  i*l,...,p  and  there  exists  a  neighborhood  VcV^  of 
throughout  which  (e^O)»  •  •  •  »®pO))  is  linearly  independent.  Restricting 
to  V,  1)  -  3)  hold. 

Since 

e^O)  £  span  [c1(uj)>...>cn_p(m)}J-, 

from  a  dimensionality  argument  it  follows  that 

span{e^(aj) , . . .  ,epO)i  =  span{c  ^  (to) , . . .  »cn_  (uo)}x  ^Ker  L(<o) 

for  all  to£v.  Finally,  let  xcX,  From  lemma  4.3  there  exist  —  C, 

i  n,  continuous  at  <o^,  such  that 

x  =  0^0^  (to)  +  . . .  +  ap(<o)epO)  +  ap+i  (|JU)C1  (u>)  +  •  •  •  +  anO)cn_p(tu) 
for  all  tu£v.  Then 


Pspan{ e L (u,) , . . . ,  ep  (cu)  } x  =  al  <“> el  ^  +  •  •  *  +  “p  eP 

30  UJ_'Pspanie1(uJ),...,e  («,)}  is  ac  “0  ^d  so  is  <o  -  spanj^U) , . . . , 

ep(uu)}.  This  completes  the  proof. 

Corollary:  Let  (v^,...,vp)  be  a  basis  of  R(uuq).  There  exist  a  neighborhood 

V  of  iu  and  maps  e.  :  Q—X,  i  =  l,... ,p,  continuous  at  tu  ,  with  e.  (u)  )  =v 
O  1  r  7  q  x  o  i. 

i*l,...,p  and  (e^(to) , . . .  ,epO))  a  basis  of  R(iu)  for  each  <ju€v. 

Proof:  Since  R  is  continuous  at  ioq,  dim  R(to)  ap  for  all  tu  in  some  neighbor¬ 
hood  of  uuq.  Setting 


LO)  =  P. 


RO)- 


and  applying  1)  -4)  give  the  desired  result. 
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We  shall  shortly  be  considering  parametrically  varying  linear 
operators.  Unfortunately , their  domains  may  also  depend  on  a  parameter. 
Hence,  in  order  to  talk  about  continuous  behavior  locally  about  the 
singular  point,  it  is  necessary  to  topologize  the  set 

Hom(R,R).  (4.2) 

To  do  this  define  :  X  -*Hom(X,X)  according  to 

I  Tx  if  x  c  R 

r (T)x  =  <  .  (4.3) 

|°  i£*'Ri 

Hom(X,X)  has  the  topology  induced  by  the  operator  norm  (3.20).  Let  i/(X) 
be  the  weakest  topology  on  X  that  makes  K  continuous.  The  topological 
space  (Z,3(X))  is  pseudometrizable  with  pseudometric 

VW  =li»-(Tl)  -►*<T2)  *  (4.4) 

Hence,  a  map  L  :  f.  —Z  is  continuous  at  if  and  only  if  ^oL  is  continuous 

at  in  the  usual  norm  sense.  Define  H(X)  to  be  the  set  of  all  L  : 

continuous  at  uj  . 

o 

If  R  :  is  continuous  at  w  we  denote  by  H  (X)  the  set  of  all 

o  R 

LSH(X)  with  L(ui)  €Hom(R(x)  ,R(uu))  for  each  u,  € If  G  and  R(ju)  *G  for 
all  uu  then  HG(X)  is  alternative  notation  for  H^(X).  Note  that  Hom(G,G)c2 
with 

iiiii  =  (it*  (T)  ii  (4.5) 

for  any  T  £Hom(G,G) .  Thus,  considering  Hom(G, G)  as  a  topological  subspace 
of  X,  the  relative  topology  on  Hom(G,G)  is  the  same  as  the  norm  topology. 
Therefore,  L  :  f.  "•  Hom(G,G)  is  contained  in  Hg(X)  if  and  only  if  it  is  con- 
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cinuous  at  <o with  respect  to  the  norm  topology  on  Hom(G,G). 

Consider  a  locally  continuous  R  :  Q  —  and  L€H  (X).  Let  V  be  a 

K. 

neighborhood  of  and  e^,...,e  basis  functions  as  described  in  the 
corollary  of  lemma  4.2.  Then 

liL(<")ei(o))  -  L(^o)ei(uuo)li  =  ||M'(L(uu))ei((ju)  -  ) e^ C^) il 


-  1;(h3L)(od)  -  (poL)  (u)o)i|i!ei(0Jo)il  +  i||i(L((Uo))ii  iie^cu)  -  e±  (ouq)  i 


(4.6) 


so  iu—  L(<o)e.(io)  is  continuous  at  <0Q,  i  =  l,...,p.  It  follows  from  lemma  4.1 
that  the  matrix  representation  of  L(w),  defined  on  V  with  respect  to 
e-.,...,e  ,  varies  continuously  with  the  parameter  about  the  singular  point. 


Proceeding  similarly,  let 


\  Hom(U.R)  (4.7) 

and  define  v  ;  l(,-Hom(U,X)  by  v(T)u=Iu  for  any  u  €u.  v  simply  extends 
range  spaces  to  X.  Letting  3(11)  be  the  weakest  topology  on  V,  that  makes 
v  continuous,  define  H(U)  to  be  the  set  of  all  L  :  Pi-’U,  continuous  at  to 

o 

with  respect  to  3(U).  If  R  :  is  continuous  at  to  ,  let  H  (U)  consist 

O  K. 

of  all  l£H(U)  with  L(to)  €  Hom(U,  R(*o) )  for  all  to.  Local  continuity  of 
matrix  representations  of  members  of  H^(U)  can  easily  be  shown. 

In  dealing  with  singularly  perturbed  systems  we  shall  consider 
only  operator  valued  maps  in  H_(X)  and  H„(U)  for  some  R,  continuous  at  to  . 
Hence,  all  systems  that  we  shall  consider  will  have  locally  continuous 
matrix  representations.  Since  matrix  representations  are  inevitably  used 
in  applications  we  are  justified  in  developing  a  theory  which  guarantees 
local  continuity  with  respect  to  the  topologies  <i(X)  and  3(U).  For 
theoretical  purposes,  however,  abstract  topological  concepts  are  preferred 
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over  notions  of  matrix  convergence  since  the  abstract  approach  is  coordinate 


free. 


We  need  a  few  more  results  concerning  H„(X)  and  HL(U).  Let 

K  R 


R  :  be  continuous  at  a  . 

o 

Proposition  4.1.  1)  Let  L  fcH  (X)  and  let  R(oj)  be  L(x)  -  invariant  for  all 

A 

xt...  Then  x— L(x)!r(x)  is  continuous  at  x  . 

o 

2)  If  X  :  U  —  C  is  continuous  at  xq  and  L &H(X)  then 

ju  —  \(x)L(x)  is  continuous  at  x  . 

o 

3)  If  £H(X)  such  that,  for  each  a)  € Ci,  L^(x)  and 

L~  (w)  share  the  same  domain,  then  tv  -  L.  (x)  +  L„  (x)  is  continuous  at  tv  . 

L  1  z  o 

4)  If  L  €H^(X)  and  L(x)  is  invertible  for  each  tvfevl 

then  x-L(x)  ^  is  continuous  at  w  . 

o 

5)  If  G  :  and  v  :  w-X  are  continuous  at  x  with 

o 

R(x)«G(x)-X  for  all  xcf.  then  M-Pj^^vW^nd  "  ~  PG(x)R(x)v(uj)  are 


continuous  at  x 


6)  If  B€HX(U)  then  x - PR(ttj)G((u)B(u>)  and  «»-Pg((b)r(u))B(u» 

are  continuous  at  x  . 

o 

Proof :  1)  Let  x£X  and  e^,...,e^  be  as  in  the  corollary  to  lemma  4.2. 

Since  pR(aj)_  3  1  “  Pj^y  -  R(x)x  is  continuous  at  x^  and  e^^ , . . . ,  e^  may  be 
constructed  to  form  a  locally  continuous  basis  for  R(x)x.  As  in  lemma 
4.1  construct  such  that 


x  =  a^(x)e^(x)  +  . . .  +an(uu)en(x) . 


^(L(x)  I  R(x)  )x  =  (x)L(x)e^  (x) 


^Pr  cHom(X,R(x))  is  the  skew  projection  operator  on  R(x)  along  G(x). 
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and  a  ~p.(L(a)  Ir('Jj))  is  continuous  at  a  . 

2)  This  follows  immediately  from 

-(X(a)L((u))  =  \(ou)^(L(oj)). 

3)  The  result  follows  from 

+L2(uj)  =  p(L^(uj))  +  t-  (L2  (ou) ) . 

4)  Construct  a  locally  continuous  na  trix  representation  with 

respect  to  a  basis  (e^,...,ep)  on  a  neighborhood  V  of  a^.  The  inverse  of 
Mat  L(a)  also  varies  continuously  at  a^.  Define  :V— C  according  to 

Mat  L(oj)  =  [3^  (a)]. 

Choose  x€x  and  construct  or  , . .  a  as  in  1).  Then 

1  P 

i.(L(a)_1)x  =  i|1ai((ju)L(uj)'1ei(u.')  =  ^(^(cb) (a))  . 

5)  Choose  basis  functions  e.,...,e  for  R  and  e  ,,...,6  for  G, 

■L  p  p*rl  cl 

and  let  a^,...,or  satisfy 


v(a)  =a1(a)e1(a)  +  ...  +  an  (a)en  (a) . 


Then 


PR(a)G(a)V(uj)  *  +  *  *  •  +  “>>ep<«>)  • 

6)  Let  u  €  U.  Then 

V(PR(a)G(a)B(Uj))u  =  PR(a)G(a)B  ('U))u- 


Letting 


v(a)  =  B(a)u  =»  v(B(a))u 


we  have  that  v  is  continuous  at  a  .  The  desired  result  follows  from  5)  and 
the  proof  is  complete. 
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4.2.  Eigenvalue  Behavior 

We  shall  now  construct  a  singularly  perturbed  system.  We  call  .1 
the  parameter  space  and  ajQ  the  singular  point.  For  applications,  the 
importance  of  non-Euclidean  parameter  spaces  is  not  as  yet  clear.  However, 
for  our  purposes,  the  structure  of  a  topological  space  will  be  sufficient. 

We  shall  not  cloud  the  relevant  issues  by  imposing  additional  structure  on 

A  generalized  singularly  perturbed  system  is  a  family  of  linear 
time-invariant  systems  described  by 

E(ou)x  =  A(oj)x  +  B(u>)u  (4.8) 

where  E,  A^H^CX),  b€Hx(U),  and  u>  ranges  over  f).  Furthermore,  we  require 
that 

det(E(<juQ)s  -  A(u>o))  ¥  0  (4.9) 

and  that  E(Wq)  be  singular. 

The  singularity  of  E(U)q)  and  the  continuity  of  E  and  A  at  wq  work 
together  to  create  the  "singular"  behavior  of  singularly  perturbed  systems. 
For  example,  we  shall  now  see  that,  in  general,  systems  of  the  form  (4.8) 
have  eigenvalues  which  can  be  separated  into  two  classes  according  to 
magnitudes  in  a  natural  way. 

Consider  the  characteristic  polynomial  det(E (uj)s-A(uu) )  of  (4.8). 
Since  forming  the  determinant  involves  only  sums  and  products  of  the  entries 
of  E(m)  and  A(uj)  we  have 

det(E(<i»)s-A(u>))  =Yn(yu)sn  +  ...  +Y1(«j)s  +  Yq(uj)  (4.10) 

where  :  Q-C  is  continuous  at  uu  ,  i*0,...,n.  Let 
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r  =maxii  I  Yi(ojo>  ^  0] .  (4.11) 

Lemma  4.3.  Lec  f  (<u,s)  =Yn(<u)sn  +  . , .  +Yq  ((jo)  where  Y  .  :  Q-*C,  i  =  0, . . .  ,n  are 

continuous  at  uj  satisfying  Y  ,  ,  )=...=  Y  (uj  )  *  0,  Y  (u>  )  t  0,  r  <n,  but 

o  r+1  o  n  o  r  o  ’ 

otherwise  arbitrary.  Then  there  exists  a  neighborhood  V  of  jjo  and  maps 
°i>  •  •  •  »an-r  ’’  V  ~‘C’  continuous  at  WQ»  with  c.(coo)*0, 
i  =  l,...,r,  ^o(ai)^0  for  all  uo,  and  such  that 

f(u),s)  ==  (uj)  C^2i  (<^i(ti))s-l))(if1(s-Xi(aj))) 


for  every  uo  €V. 

Proof;  Let  g(uo,s)  =  Y0(m)sn  +  . . .  +  Yn(uo) .  A  bound  on  the  roots  of  a  polyno 
mial  over  C  given  in  [26],  p.  62  implies  that  there  exists  at  least  one 
root  6^  of  g (uu, s )  satisfying 


*Y  (a») 


iJi<^¥) 


i 

a-r 


(4.12) 


whenever  Yr(<u)  ^  0.  Since  Yr(a)Q)  ^0,  there  is  a  neighborhood  of 
throughout  which  (4.12)  holds.  Define  on  according  to 


aL(ui)  =*  / 


6 

OJ 


if  Yn(cu)  *0 


^0  ifYn(aj)=0 

The  construction  of  guarantees  that  if  0  is  a  root  of  g(uo,s)  then 

a,  (<u)  »0.  From  (4.12)  a,  is  continuous  at  oo  with  a,  (uo  )  =0. 

1  1  O  i  o 

Let  C(VL)  be  the  set  of  all  maps  from  into  C,  continuous  at 
uoq.  Using  pointwise  addition  and  multiplication,  C(V^)  is  a  commutative 
ring  with  identity.  Hence  (see  [25], p.  334)  there  exist  YQ, . • • >Yn_^  6C(V^) 


such  that 
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g  =  (s-a1)(Yosn'1  +  ...+Yn.L) 

where  g  is  considered  as  a  polynomial  over  C(V^).  Multiplying  and  equating 

coefficients  at  ou  yields  Y  .  (cu  )*...*Y„  ,  (uj  )  =  0  and  Y  (uj  )  t  0. 

o  r-t-i  o  n-  j.  o  r  o 

The  above  arguments  may  be  applied  n-r  times  yielding  a  neighbor¬ 
hood  V  of  tu  and  maps  o  , ...  ,a  ,  0  , . . .  ,<p  :  V  — C,  continuous  at  uu 
n-r  o  l  n-r  or  o 

with  c.  0"  )*...=  3  (cu  )  =  0  such  that 
L  o  n-r  o 

g=  (-l)n_r(^rsr  +  ...  +^0)^(s-ai)- 

From  the  construction  of  the  a.  we  have  that  0  is  not  a  root  of 

x 

2^(00)  sr  +...+  g>Q  (uu)  so  $q(uu)#0  for  any  ,JU  ^vn_r*  Let  v  =vn  r’ 

Clearly, 

f(uu,s)  =  (vpQ (jo) sr  +  . . .  +  *>r(JJ^i=1_(<Ji(tJ,)s"l)  * 

From  the  continuity  of  a  and  it  follows  that  there  exist  :V-*C, 

continuous  at  uu  such  that 
o 

r  r 

^o(uu)s  +  . .  .  +0r(u))  =  ^Q(<u)i2]_(s  '  ^0*0) 

for  all  uu£v.  This  completes  the  proof. 

We  thus  have  that  the  eigenvalues  of  the  system  (4.8)  are 


X1(ui)J...»X.r(UJ)  atld  a  _  (oj)  ’  ‘  • '  ‘  a  (w) 

1  n-r 


where  and  a.  are  continuous  at 


0)^  and  a^(0UQ)  =  0.  Hence  there  is  a  characteristic  separation  of  the  modes 
of  all  systems  of  the  form  (4.8)  according  to  magnitudes. 

It  is  now  a  straightforward  task  to  prove  proposition  3.1. 

Proof  of  Proposition  3.1.  Let  0  =  ?^  with  the  topology  induced  by  the  norm 
(3.20)  and  choose  $  J.  Then  from  lemma  4.3 

.  P  n'r‘P 

det(Lfs  -  I  -  Bf(»  I  F))  =  «(«)  (.^  (s-,,.  (uu) ))  (^  (6.  (co)s-l)) 
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for  some  p£n-r  where  a(u)Q)  ^0  and  6^(00)  =0.  Clearly, 

det(Lfs  -  I  -  Bf  (a)  I  F))  ¥  0 
is  satisfied  throughout  a  neighborhood  of  ujq.  Also 

£a(E,A) 

implies 

\(»>  t  a(E,A) 

throughout  a  neighborhood  of  ujq.  Finally, 


1 

6i((ju) 


jE  ^(E,A) 


so 

tj(Lf,I  +  Bf(oa  |  F))  fi  cr(E,A)  =0 

throughout  some  neighborhood  of  u)q.  Hence  T  is  open  in  Q  and  the  proof  is 
complete. 


4.3.  Slow  and  Fast  Subspaces 

Since  the  system  (4.8)  is  a  descriptor  system  at  U)q,  the  slow 

and  fast  gubspaces  are  already  defined  at  m  .  We  shall  now  extend  their 

definitions  to  a  neighborhood  of  u>  .  From  lemma  4.3  there  exists  a 

neighborhood  V  of  m  such  that  | X .  (tu)  |  <| — 7—rl  for  all  uu€V.  We  now  re- 
or 

define  Q  =  V  so  that  the  natural  eigenvalue  separation  occurs  at  all  points 
of  G.  Since  we  are  concerned  only  with  local  properties  of  (4.8)  about 
the  restriction  of  Q  to  a  neighborhood  of  <juq  causes  no  problems. 

To  facilitate  the  definition  of  slow  and  fast  subspaces  we  need 


the  following  lemma 
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Lemma  4.4.  Let  E,  ASH^(X)  satisfy  (4.9).  Then  there  exists  X 
continuous  at  id  such  that  X(x)  €c  -  a(E(ui),A(oj))  for  all  u) 

Proof:  Let  V,  :fi— C  be  continuous  at  oj^  with  ri(^0)  it  c(E  (u)q)  ,A(uuo) ) .  Then 
there  exists  a  neighborhood  V  of  id  such  that  ‘ ;  (oo)  ^  X^  (oj)  ,  i  =  l,,..,r  and 
1  V|(aJ)  |  <  | -y— | ,  i  =  1, . . .  ,n-r.  Hence  "(ui)  £  a(E  (ui)  ,A(>d) )  foralliutV. 
Define  X («JJ)  =  .,(«>)  for  oj  €v  and  let  X(w)  satisfy  X(w)£a(E(j), 

A(x) )  for  <u£ft-V.  This  completes  the  proof. 

Choosing  X  as  in  lemma  4.4  we  may  now  define  the  slow  subspace 

at  oi  as 


S(w)  =Ker  (X  (uu)E  (m)  -A(cu))_  E(u>)  -■ 


and  the  fast  subspace  at  o)  as 


.(x)  -  Xi(tu) 


cr .  (x) 


(4.13) 


F(»)  =Ker  ^((X  («)£(«)  -  A(cd)  )-LE  (ao)  -  — 1) .  (4.14) 


1  r  JL)  . 

det (Is  -  (X(x)E(x)-A(W))'LE(u)))  -^(s  *  \J  1 


(4.15) 


with  i  ^  j  implying  ^T.  .  As  demonstrated  in  the  proof  of  theorem  2.1, 

im  joi  r  5 

X (uj)  -X  (oj)  an  eiSenvalue  (X(<u)E(u))-A(«u))  ^E(oi)  and  so  may  be  identi¬ 
fied,  after  proper  indexing  with  ™.  ,  i*l,...,q  for  some  q  .  The  re- 

1<JJ  .1  ai^) 

maining  eigenvalues  of  (X(o))e(<d)-A(id))  l(ui)  are  ^  and  may  be 

identified  with  7)^,  i  = q^+l, . . . , p^.  Hence  we  may  give  an  alternate 
definition  of  S(<d)  and  F (o))  as  the  direct  sum  of  eigenspaces.  Note  the 
resemblance  to  (2.7)  and  (2.8), 

n , 

S(u>)  =■  if1  Ker  ( (X (oi)E (x) -A (oi) ) "  t  (x)  -  H^I)  lu;  (4.16) 

Po  .  n . 

F(»)  =  i4  -  "  I)  LJ)  (4.17) 


1 
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Clearly,  the  definitions  (4.16)  and  (4.17)  are  equivalent  to 

(4.14)  and  (4.15).  (4.14)  and  (4.15)  are  usually  to  be  preferred,  however, 

since  they  are  given  in  terms  of  operators  that  vary  continuously  with  oi 

at  v)  .  It  follows  immediately  from  lemma  2.1  that  S(m)  and  F(ou)  are 
o 

independent  of  the  particular  choice  of  the  function  X  for  all  ui  60.  From 
the  definitions  and  the  properties  of  X^  and  as  outlined  in  lemma  4.3 
it  is  clear  that  dim  S(m)  =r  and  dim  F(m)  =n-r  for  all  oi.  Also,  S(iu) -®F(uj) 
=  X  for  all  tu.  From  (4.13),  (4.14),  and  lemma  4.2  we  have  that  S,  F  :  Cl -771 
are  continuous  at  o>  . 


4.4.  System  Decomposition 

We  now  extend  the  decomposition  for  descriptor  systems  to  a 
neighborhood  of  the  singular  point.  We  first  extend  the  algorithm  (2.10)  - 
(2.13)  to  all  of  Q.  Let  J^,  J2  :  Cl-Z  be  defined  by 

J^w)  =  (X(tu)E(u))-A(uu))  E(«)  I  S(u>)  (4.18) 

J2(uj)  =  (\ (uj)E (<u) -A (m)  )  _1E (w)  I  F(ou) .  (4.19) 


Define  M,  M  :  Q-*Hom(X,X)  by 

( 


-1 


_  I J.  (m)  x  if  x  6  S  (u>) 

M(oi)x  =  ly 

(\(uJ)J2(m)-I.)_1x  if  x€F(<u) 


and 


(4.20) 


M(cu)  =M(a>)  (X(cu)E  (m)-A(cu)) 


Finally,  let  L  ,  :  Q-Z  be  given  by 


(4.21) 
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Ls(oj)  =\(uj)I  -  (4.22) 

and 

Lf(<u)  =  (X(aa)  J2  (oj)-I>'LJ2  (uj)  .  (4.23) 

Clearly,  M(o.)  and  M(x)  are  invertible  for  each  uu  6  0.  From  proposition  4.1 
we  have  LgcHs(X),  J2>  Lf  c  Hp  (X) ,  and  M,  M€HX(X). 

At  this  point  we  are  ready  to  state  the  main  decomposition  result. 
Theorem  4.1.  For  each  uu 

1)  S(a>)  and  F (uu)  are  both  M(uu)E(uu)-  and  M(uu)A(u>)-invariant. 

2)  M(uj)E(<u)  I  S(uu)  =1,  M(uu)A(uu)  I  F(uu)  =1. 

3)  M(oj)E(uu)  I  F (uu)  =  Lf  (uu)  ,  M(cu)A(uu)  1  S  (uu)  -  L  (<fl)  • 

4)  det(Lf  (uj)s-I)  =  i^(ct(uu)s-l) . 

5)  det(Is-Lg(<JJ))  =  i21(s-Xi((u)) . 

Proof:  1)  -  3)  follow  immediately  from  the  definitions  and  from 

(X(tu)E(uu)  -A(uu))-1A(uj)  =  \(uu)  (X(uu)E(’ju)  -  A  (uu)  )  ‘  h  (u))  -  I. 


.To  prove  4)  note  that 


n-r 


det(ls  -  J2(uu))  -\(UJ)ai(uu)  -  1  )- 

Then  from  (4.23)  we  have  the  desired  result.  5)  follows  similarly  by 
observing  that 


det(Is  -  V®»  "  iVs  -I^)TW  >• 

This  completes  the  proof. 

From  lemma  2.2,  M  is  independent  of  \  and  hence  so  are  L  and  Le. 

S  II 

From  6)  of  proposition  4.1  we  may  define  B„  €HC(U)  and  Bp  €H  (U)  according 

s  b  I  r 

Co 
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Bs<“>  =  PS(.)F(.)BW  <4.24) 

and 

sf(“>  ■  'iwswX’*'  <4-25> 

We  have  finally  arrived  at  the  decoupled  system  equivalent  to  (4.8).  Letting 
the  solution  x  of  (4.8)  be  decomposed  by  S(u>)  and  F(<u)  into  x=xs+x^  we 
define  the  slow  subsystem 

*9  -L8(u>)xg+Bs(u))u  (4.26) 

and  the  fast  subsystem 


(uj)x^  =  4-  (ui)u . 


(4.27) 


The  initial  condition  xq  is  decomposed  into 


x  (ui)  =  P„  ,  .  _  .  ,x 
so'  S(u>)F(u>)  o 


(4.28) 


and 


!tfo<‘",'pF(w)S<«)V 


(4.29) 


From  5)  of  proposition  4.1,  xgQ  and  x^q  are  continuous  at 

We  have  thus  achieved  an  exact  decoupling  of  (4.8)  according  to 
eigenvalue  magnitudes.  Continuity  at  the  singular  point  has  been  preserved 
everywhere  possible.  The  decomposition  at  coincides  with  the  descriptor 
decomposition  developed  in  Chapter  2.  We  shall  use  the  decoupling  results 
in  a  variety  of  situations  in  later  chapters. 
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CHAPTER  5 

TRAJECTORY  CONVERGENCE 

5.1.  Solution  of  Descriptor  Equations 

In  this  chapter  we  shall  consider  two  questions.  The  first 
concerns  one  justification  for  calling  (2.24)  the  solution  of  (2.20).  Con¬ 
sider  the  class  of  all  singularly  perturbed  systems  whose  system  operators 
converge  to  those  of  a  given  descriptor  system.  Given  one  member  of  that 
class,  its  solution  may  or  may  not  converge  as  We  would  like  to 

find  the  set  of  all  possible  limiting  solutions  of  such  systems.  If  there 
exists  only  one  possible  limit  then  it  would  make  sense  to  call  that  the 
solution  of  the  descriptor  equation.  Since  each  descriptor  system  is  the 
limit  of  a  singularly  perturbed  system,  it  would  be  convenient  if  at  least 
the  possibility  of  convergence  of  solutions  existed.  Of  course  we  have  to 
decide  what  definition  of  convergence  of  functions  to  use.  It  would  be 
helpful  if  more  than  one  notion  of  convergence  gate  the  same  result.  ‘  So  far 
no  other  author  has  taken  such  an  approach  to  the  solution  of  (2.20). 

The  second  problem  is  that  of  determining  when  the  solution  of  a 
singularly  perturbed  system  does  converge.  More  will  be  said  about  this 
later. 

We  first  consider  generalized  functions  (see  [24])  from  [0,°°) 
into  X  and  U.  Let  K(X)  be  the  C-vector  space  of  infinitely  differentiable 
maps  from  (-00,00)  into  X  having  compact  support.  A  sequence  ($,  )  in  K(X) 

K 

converges  to  j€K(X)  if  and  only  if  the  supports  of  all  the  are  contained 
in  the  same  bounded  interval  and  ^  uniformly  for  i  =  l,2,3,....  The 
class  of  generalized  functions  is  the  set  of  all  continuous  linear  func¬ 
tionals  on  K(X).  We  shall  consider  only  a  subspace  &o(X)  of  the  generalized 
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functions.  *o0O  consists  of  all  x  with  (x,£)  =0  whenever  0  h.t.s  support 
contained  in  (-^,0],  In  other  words,  ^q(X)  consists  of  all  x  such  that 
x  =  0  on  0) . 

If  x  :  [0,x)-X  is  Lebesgue  measurable  and  integrable  on  any 
compact  interval  then 


(x,£) 


•  CO 

W  o 


defines  the  corresponding  xSi^CX), 
by 


<x(t),0(t)  >dt  (5.1) 

For  v  €X,  the  functional  6Xv  defined 


(6Lv,0)  =  (-l)1<v,/(0)  > 


(5.2) 


is  the  multivariable  generalization  of  the  Dirac  delta  (differentiated  i 
times).  We  shall  consider  the  two  most  common  topologies  on  ^q(X),  the 
weak  and  strong  topologies  (see  [24]).  A  sequence  (x^)  in  **0(X)  is  said  to 
converge  weakly  to  x€4q(X)  if  (x^,sS)  -  (x,^)  for  every  £6K(X).  We  need 
not  even  consider  the  definition  of  strong  convergence.  For  our  purposes 
it  suffices  to  state  that  strong  convergence  implies  weak  convergence  and 
the  two  limits  are  equal^. 

Consider  the  descriptor  variable  system  (2.5)  and  the  class  of 
all  singularly  perturbed  systems  converging  to  it  in  the  sense  of  its  de¬ 
fining  operators.  We  could  attempt  to  show  that  there  is  only  one  limit 
for  the  solution  to  achieve,  but  we  are  already  satisfied  with  ,ie  defini¬ 
tion  of  the  solution  of  the  slow  subsystem  (since  this  is  merely  a  state 
equation).  In  fact,  only  the  natural  response  of  the  fast  subsystem  is  in 
question  since  the  forced  response  is  the  solution  of  (2.20)  in  the  ordinary 
sense  for  u£C^(U).  We  especially  need  to  justify  our  use  of  the  unforced 

^The  three  topologies  on  K(X)  and  J»o(X)  are  all  Hausdorff  and  satisfy  the 

first  axiom  of  countability.  Hence,  we  need  consider  only  countable 

sequences.  Also,  no  sequence  has  more  than  ona  limit. 
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part  of  (2.24)  since  it  does  not  even  satisfy  (2.20). 

Let  and  (T^)  be  in  £  with  invertible  in  Hon^R^R^ 

and  -  T  cHom(R,  R) ,  T  nilpotent  with  index  q.  Let  v^-v  £R.  We  need  to 
consider  the  solution  of 


T,  x  =  x 
k 


(5.3) 


with  initial  condition  v.  . 

k 

Theorem  5.1.  If  e(Tk  Svk  converges  weakly  to  some  limit  then  that  limit 

\\i-l  i 
is  T  v. 

Proof;  We  have  for  €K(X) 

-1  i’03  ^k  ^  i  i  i-L 

(e(Tk  )vk,sZ>)  =jQ<e  K  vk>0(t)  >dt  =  .^(-l)1  <1^,/  "(0)  > 


CT, 


-1 


+  (-l)qJ”<Tje  k  vk,?>q(t)  >dt 


by  integration  by  parts.  Choose  an  orthonormal  basis  (e^,...,en)  of  X  and 


let 


CTk 

6  vk  ~  alk(t)el+  +0nk(C)en 

^(t)  =  ^1(t)eL+  ...  +0n(t)en 


^(Tk  )ej  ~  ~lijkel+ •*  *  + 


bjk  n’ 


Since 


^(\q)  =  (^(Tk))q» 

'(T,q)  “*0  and  T,.  .  —  0  as  k—  ®.  Since  pq  6K(X),  we  may  define  .  €K(X)  by 
K  1J 


f .  . (t)  =  P.  (t)e  . . 

ij  v  '  L  J 


Then 
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tT'L 

;Vk(t)3(t)dt=;“<e  k  vk^ij(t)>-c.. 


for  some  C..,  i,j=l,...,n  and 
«;*J 


tT 


-1 


<Tke  k  >dt  =iiLJo<ri(Tkq)(e  K  V’3i(t)ei  >dt 

^  '.03 

=  .  5  -7).  a..  (t)£.  (t)dt  —  0. 

i, j  =  l  Ljk"  o  jk  "iv 


Hence, 


q-1 


q'La-i„i 


(e(T“1)vk,0)  -  V(- 1)1  v^v.^'^O)  >  =  (-^■Yv.W 


for  all  0  6K(X)  and  the  proof  is  complete. 

We  therefore  know  that  the  only  possible  limit  of  the  solution 
of  the  fast  subsystem  (4.27)  (unforced)  in  the  weak  topology  on  &q(X)  is 
in  fact  what  we  have  been  calling  the  solution  all  along.  To  reinforce 
this  convergence  argument,  consider  the  strong  topology  on  &q(X).  If  the 
unforced  solution  of  (4.27)  converges  in  the  strong  sense  then  it  must 
converge  weakly.  Since  the  two  limits  must  be  the  same.  Theorem  5.1  holds 
for  strong  convergence  as  well.  The  labeling  of  (2.24)  as  the  "solution" 
of  the  fast  subsystem  is  inescapable.  It  should  be  stressed,  however,  that 
(2.24)  can  be  called  the  solution  only  in  the  limiting  sense.  Once  again, 
it  does  not  satisfy  the  equation  (2.20). 


5.2.  Sufficient  Conditions  for  Convergence 

Although  the  previous  section  shows  that  only  one  limiting  solu¬ 
tion  of  a  singularly  perturbed  system  can  exist,  the  question  of  whether 
or  not  that  limit  is  actually  achieved  is  still  unanswered.  In  this  sec¬ 
tion  we  shall  develop  conditions  which  guarantee  convergence  of  the  solu¬ 
tion  as  ou  —  jj  . 

o 
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For  several  reasons  it  is  important  to  know  if  the  system  trajectory 
converges  with  tu.  Since  in  practice  a  descriptor  system  can  usually  be 
viewed  as  the  limit  of  a  particular  singularly  perturbed  system,  establish¬ 
ing  trajectory  convergence  serves  as  justification  for  considering  a  des¬ 
criptor  system  as  a  viable  model.  A.  descriptor  system  which  has  not  been 
identified  as  the  limit  of  a  particular  singularly  perturbed  system,  or  for 
which  trajectory  convergence  has  not  been  established  may  have  hidded 
instability  due  to  the  disappearance  of  the  fast  modes  in  the  limit. 

In  terms  of  numerical  considerations,  it  is  easier  in  some  cases 
to  compute  the  solution  of  an  rth  order  descriptor  system  than  an  nth  order 
state  variable  system.  If  trajectory  convergence  has  been  established,  a 
computational  savings  may  be  obtained  by  setting  the  parameter  to  tu  and 
calculating  the  corresponding  solution,  sacrificing  a  small  amount  of 
accuracy. 

Convergence  of  the  solution  of  the  standard  system  (1.5)  is 
understood  (see  [20], [28],  and  [34]-[36]).  Necessary  and  sufficient  condi¬ 
tions  for  existence  of  a  limiting  solution  can  be  stated  in  terms  of  the 
limiting  system  alone.  This  can  be  done,  however,  only  because  the  assump¬ 
tion  of  the  form  (1.5)  contains  implicitly  the  assumption  of  one  particular 
way  of  approaching  the  limiting  system.  In  our  formulation  many  different 
approaches  are  possible  leading  to  a  much  more  difficult  problem.  To 
guarantee  convergence  of  solutions,  some  statements  must  be  made  not  only 
about  the  limiting  system,  but  also  the  way  in  which  it  is  approached. 

Convergence  of  the  solution  of  the  slow  subsystem  (4,26)  in 
various  senses  is  relatively  easy  to  establish.  We  shall  be  concerned 
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mainly  with  characterizing  convergence  of  the  solution  of  the  fast  sub- 

1 


system  (4.27)  in  terms  of  the  behavior  of  the  fast  eigenvalues 


a.(uj)  * 


For  convenience  we  assume  that  Lc(u))  is  invertible  when  .  That  is, 

r  o 


is  an  isolated  singularity. 


Let  $£  :  Q  x  [O,00)  x  c"(U)  xX  —X  be  given  by 


/  tLf (uj)  ,t  (t-T)Lf(a))'L 

e  xfo(u))+Je  Bf(u))u(T>dT  if  uj^iu 


§f  (ai,t,u,xo)  =  ( 


q-i 


(5.4) 


-.2  L  (in  )1B,(uj  )uX(t)  if  'ju=oj 
^  1=0  fv  o  f  o  '  o 

where  u1(0)  is  the  ith  right-hand  derivative  of  u  at  0.  Clearly,  $ 
agrees  on  (0,°°)  with  the  solution  of  the  fast  part  of  (4.8)  at  uj  with  con¬ 
trol  u  and  initial  condition  x  . 

o 


From  (5.4)  it  is  clear  that  understanding  the  behavior  of 


-1 


e(L^(yu)  )  plays  an  essential  role  in  the  study  of  convergence.  We  now 
characterize  convergence  of  e(L^(tu)  )  in  terms  of  the  fast  eigenvalues 


o.  (<jj)  ' 


Consider  the  smallest  rectangle  in  C,  symmetric  about  the  real 


axis,  enclosing  the  eigenvalues  "  ■.  y  .  (See  Figure  5.1.)  Let 

j  ^  ' 


a(ui)  =  max  Re 


b(«)  =  min  Re 


°j  («) 

1 

cr.(o)) 


c(uj)  =max  |  Im 


(<U) 


(5.5) 

(5.6) 

(5.7) 


Define 


“Here  "Im"  denotes  imaginary  part. 
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Y  (w)  =  mini. 


Ia(a»> 


2  2 
a('JU)  +c(cu) 


bWl  } 

)  1  J  • 


b  (tu)  ^  +  c  (tu)^ 


(5.8) 


As  a  part  of  any  sufficient  condition  for  convergence  that  we 
shall  derive,  it  will  always  be  assumed  that  a(<u)  —  as  UJ-‘U,0  and  hence 
b(<u)  —  -00.  Actually,  it  can  be  shown  that  this  is  in  fact  necessary  for 
most  types  of  convergence.  Thus  we  may  restrict  attention  to  a  neighborhood 
of  cuq  with  iu€v^  -  {(Uq}  implying  a(to)  <0.  Define  p^  :  [0,1]  ~*C  by 


P  (y)  a  — |ei2TTy  .  (X££2+_i — ) 

Pujly;  i  2  a(tu),e  ^  Y  aOu)'* 


(<u)' 


(5.9) 


p^  parametrizes  a  circle  with  center  at 


and  radius  |M+  1 


a(uj) 


a  (<u) 


Let  p  =6  +V|  i. 
r'JU  O)  ou 


Lemma  5.1.  o^y)  <0*  aj  is  encl-osed  by  P^j  and 

Ip  (y)  -cri(UJ)! 

0)  J  I 


for  j  =  l, . . .  ,n-r  and  each  y  €  [0,1],  to  €  V ^  -  [u)q}  . 
Proof :  Choose  j,y,  and  tu.  We  have 


From  (5.9), 


Y(uj)  * 


-a(m) 


a(uj)^  +  c  (iu)^ 


_1 _ 

a(uj) 


<  0. 


Let  a^(rn)  =w  +  zi. 


Then 


bWiy-j  *a(<») 
v  +  z 

and 

£  c"")- 

w  +  z 
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Equivalently, 


12  2 
(W-2lW}  +Z  * 


1  ,2  2 

(w-2b^)>  +2  " 


2  12 
w  +  (z  -  r  "7  --)  2: 
v  2c  (uu)' 


2  12 

"  +(s  +  2^>  2 


1 

4a(uj)2 

1 

4b  (uo)2 
1 

4c  (ou) 2 
1 

4c  (uo)2 


(5.10) 


Hence,  a.(ui)  is  contained  in  the  region  determined  by  four  circles  as  in 

c  i  -m.  c.  .  ^  .  .  a  (oo)  +  c  (tu)  i  b  (uo)  ±  c  (uQi  „ 

Figure  5. The  five  points  of  intersection  are  ■■  ~ — N  ,  — ^2 -  — 2’ u* 

a(uo/  +  c(u>)2  b(u))  +  c(w) 

From  the  geometry  it  is  clear  that  (ou)  .  From  (5.10)  it  follows  that 


w  +  z  s 


w 


a  (<*>) 


so 


(w 


.Yiii._±_)2  +  Z2£  +maii+-± 

4  a(uu)  '  2  a<',n'v;W+  4  +af-.;V  “  16  +  ?a  r  + 


a(u>) 

=  (3IM+JL_)  2 

^  4  a  (uu)  ^ 


a(ti)- 


16 


2a(u)  ,  ,2 

a(oi) 


Thus,  a\(uo)  is  enclosed  by  the  circle  with  center  at  -^7-^  +  7  -  and  radius 
J  4  a(uu) 


a(yu) 


a(uj) 


and  the  proof  is  complete. 

Lemma  5.2.  If  a(uo)--®,  b  (uu)$a  ^  -  0,  and  c^)^3^  -0  as  uu-cu  for 
-  0 

every  p  >  1  then  e(L^(uo)  S~0  pointwise  on  (0,33)  and  uniformly  on  [€,=°) 
with  respect  to  the  pseudometric  d  (see  (4.4))  for  all  €>0. 


FP-6774 


ntaining  ck  (uj)  . 


\ 


r 


Proof:  Let  tsi>0,  Then 


,  -1  t 

CL^(0))  J.  .  g  -l 

e  =‘j~r  $>  e  (sI-Lf(x))  ds 


where  the  path  of  integration  is  parameterized  by  p  .  Hence, 

t  ^ 


tLf(x) 


P.Jy)  p  ’  (y) 

L  ,  1  e  x 


"  2^'  J  o  n-r 


"  (ad j  (p  (y)I  -  Lf  (aj)))dy  : 

dl  £ 


.1,  (p  (y)  -  -  ,(■-«) 

J  — L  (JO  J 


-L|ePa,(y)|«| K£h 

“2TT  o  Wp 


jSi'p^y)  -  V^1 


-  f*(adj(p  (y)I  -  Lf(uu))):!dy 

'L  £ 


where  "adj"  denotes  the  (classical)  adjoint.  Since  p  (y)— 0  as  uj-uj^ 


uniformly  in  y, 


f* (adj  (p  (y)I  -  Lf  (oi) ) )  —  -t~ (adj  Lf  (uj  ) ) 

UJ  £  £0 


uniformly  in  y  and  there  exists  a  neighborhood  V  of  cjuq  and  N>0  such  that 


ladj  (P^Cy)!  -  Lf  (uj)  ;!  <N 


for  all  uj  € V  -  {uj  } ,  y  t [0, 1] 
We  have  for  uj  £  -  Iuj^}  , 


(-<Vy;  4  a  (uj)  '  ^  V7'  Kl  a(u;);  a(uj); 


Vy)  1  „  2  .  . 

.  ,  N2  ~  .  ,2  -  .  ,Y(uj)  ^  _1  .  "  3 
5  (y)  +  ~l  (y)  2  (  f  ■  *  +  -7— r) 

uj  7  ujw'  v  4  a(uj) 


since  Y (uj)  S  - — 7— r.  Thus 
a(uj) 
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l.,.<y,l..3€*w1 


From 


ip;(y)|.2n|iM  +  ^| 


it  follows  that  for  id  €v  Hv.  -  Icd  J 

1  o 


-1 


2  . 


j|KeCLf(U0 


n-r 


which  is  independent  of  t.  If  n-r  =  1  we  need  only  show  that 


Letting  $  =  e  we  have 


j  ca(«D) 

e  _ 

a(u))V  ((d) 


2a  (id) 


a  (id)  ( - - 5) 

afu)-)  '  +  c((D)“ 


2a ((D)  .1  ,  ,  v  ^a((D)  „ 2 

=  -P  + - r  (c(cd)0  v  ) 

a((D)z 


and 


2a  ((D) 


a  ((D)  ( - - t)  a  (id) 

b(<D)  +c(<d) 


1  j((b((D)«a(<U))2  +  (c(<D)0a(tU>)2). 


For 


n-r  >  1  set  0  =  e"^n  r  ^  .  Then 


6  a  ((d)  2a((U)  n-r-1 

6  "  (^™-) 

Y((D) 


n-r-1  v  V  ((D) 


But 


2a  ((D) 


2  2 
a  (<d)  +  c  ((D) 


a(uu)02a(U))  +-A^(c((D)0aW)2 


75 


*2a(») 

b(m) 

2  2 
b(ou)  +c(ou) 

€ 


b^y((b(o))/(u,))2  +  (c(<u)/(lu))2). 


Finally,  setting  £=e2^n  and  applying  similar  arguments  gives 


|€a((U) 


a(ua)Y  (<o) 


This  completes  the  proof. 

Lemma  5.2  establishes  convergence  of  the  exponential  provided 
a(o>)  —  -00  and  certain  growth  conditions  on  b(iu)  and  c(ou)  are  satisfied.  For 


example,  if 


b(u>)  a  (id) 


(5.11) 


for  some  N,  k>0  then  the  condition  b(v))0a^aJ^  — 0  holds  for  every  <t>  >  1.  The 
rate  of  growth  of  b(<ju)  and  c(<d)  must  be  less  than  any  exponential  of  a(uj). 
Lemma_5JL3.  If  a(OJ)--00  then  there  exists  a  neighborhood  V  of  such  that 


co  ,  tLf(uj) 

i  t 


!|dt  < 


Ib(uu)2  +C(U))2 


for  all  ou  £V  -  U  } . 

o 

Proof :  As  in  the  proof  of  lemma  5.2, 

_1  — - — 

tLc(m)  .  n-r  v .  ,  1  o  (y) 

Me  f  ||  *»(-*-)  l^  +  -M  i  le®  |dy 

1  11  ',Y(UJ);  2  a(u>)  q  1  y 

so,  by  Fubini's  theorem, 

tLf(u,)~\  4  n-r  |Y  ((d)  1  |  1  °?  P^, 

0  :ie  !'dtS:N(Y^))  ,“P+Iwl/o;Je  |dC  dy 

O 

£fj(- — — — )  (T.W  + — i — w 2 )  S  — N  ( — — — )  ( — i — )2 

'  2  a(<D)M2a(uV  2  va  (ou)  ^ 
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Since  a (oj)-'-03,  there  exists  a  neighborhood  V  of  ou  with 

( — " — <  —  (■ — - — ) 

ka(uJ);  3  ,  n-r.. 


Then 


4“  N 


=o  tL,(tuf^  ,  n-r  ,  ,  ,2  ,  .2 

,  !le  f  llHt.<  /_!_>  gib(qj)  +c(o>)_ 

Jq"  dt  S' (cur  a('ju) 


n-r 


and  the  proof  is  complete. 

We  now  establish  sufficient  conditions  for  convergence  of 

Theorem  5.2.  Let  a(U))--x,  b  (cu)0a  ^  -  0,  and  c(<u)0a^-O  for  every  <t>  >  1 
00 

and  let  u  $C  (U)  have  all  its  derivatives  bounded.  If  there  exists  a 
neighborhood  V  of  <u^  and  a  positive  integer  p  such  that 


!|L,(uO  ii  < 


for  all  uu  £v  -  U)  } ,  then,  as  <u  -<v  , 
o  o 


Jli 


n-r 


b(W)2  +  c  (ul>)  2 


«f(0), -,u,xo)  -4f(u,0>*»u,xo) 

pointwise  on  (0,=°)  and  uniformly  on  [£,=°)  for  every 

c>0  and  each  x  £x. 

o 

Proof:  For  cu^cu^,  integration  by  parts  gives 

tLf((u)_1  P  i  i 

}f  ((U,t,u,xo)  =  e  xfo(ou)  -  i^0  Lf  (<u)  Bf(<«)u  (t) 

tLf(w)  1  P  ±  ±  ^(t-t  )Lf  (d»)  1 

+  e  1|0Lf(tu)  Bf(o))u  (0)+Lf(<u)pje,  Bf  (cu)uP(T)dT. 


From 


aW 


.  ,  v2  ,2'  'b(uu) 

b  (cu)  +  c  (cu)  v 
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p 

it  follows  that  Lf(w)  -0  so  paq.  Hence,  from  lemmas 
desired  result  follows. 
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CHAPTER  6 

FURTHER  STRUCTURAL  PROPERTIES 


6.1.  Slow  Subsystem 

It  is  often  advantageous  to  know  that  if  a  certain  dynamical 

system  property  holds  for  a  singularly  perturbed  system  at  m  ,  then  it  also 

■  '  o 

holds  throughout  some  neighborhood  of  u)  .  For  example,  in  this  chapter  we 

shall  see  that  controllability  behaves  in  this  way.  Hence,  controllability 

at  is  sufficient  to  guarantee  existence  of  the  optimal  solution  of  the  LQ 

regulator  problem  (see  Chapter  7)  throughout  a  neighborhood  of  uu  .  This 

o 

fact  will  result  in  a  computational  savings. 

In  this  chapter  we  shall  consider  various  results  of  tnis  type. 
Applications  to  the  pole  placement  problem  will  be  discussed. 

In  section  4.3  we  restricted  the  parameter  space  0  to  a  neighbor¬ 


hood  of  ai  so  that  the  eigenvalues  of  the  system  4.  (u))  and  are 

O  1 

separated  by  magnitudes  for  all  uj  fcQ.  Since  each  X^  is  continuous  at 

the  parameter  space  can  be  further  restricted  so  that  ^(u)Q)  implies 

X^(<jo)  A-j  (tu)  for  all  uu  6Q.  Accordingly,  the  functions  V  can  be  partitioned 

into  equivalence  classes,  X^  and  X^.  being  in  the  same  class  if  and  only  if 

(<ju  )*X.(to  ).  We  can  reindex  the  X.  such  that  X.,...,X  comprise  one 

iv  o  j  o  i  1  n^  r 

equivalence  class,  X  X  another,  etc.  up  to  n,  =  r. 

Consider  the  eigenspace 

sl<“>  ■Keti-n"1l.l+l<a<'")El“>'A<”))'LE(“>  '  IWXW  I} 


n . 


(6.1) 
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} 


where  a  =0  and  \  is  as  in  lemma  4.4.  From  lemma  4.2  and  the  fact  that 
o 

S^(uj)  has  constant  dimension  throughout  fl  it  follows  that  is  continuous 
at  for  i  =  l,...,k.  Also, 

k 

SO)  =  ,f1Si(a»).  (6.2) 


Lemma  6.1.  Let  R-.....R  :  Q-'77(  be  continuous  at  oi  .  Then  there  exists 

-  ]_»  p  o 

R  ;  Q  “*>'/(  such  that 

1)  R  is  continuous  at  ou  . 

o 

2)  R(uu)CR1(uj)  +  ...  +Rp(tu)  Vo)€Q. 

3)  R(uio)  -R1Oo)  +  ...  +Rp(u»0). 

Proof ;  Let  L  :  Q  - Hom(XP,X)  be  defined  by 


LO)(x1,...,Xp)  ■PRi(0U)*1+  •••  +pRp(u,)Xp" 

Then  uu-*L(<JJ)  is  continuous  at  and  by  lemma  4.2  there  exists  R  :  Q  —  //J, 
continuous  at  <jjq,  with  R(uu)  ^Ker  L(uj)*  for  all  uj  €0  and  R(u)q)  =  Ker  LOo)*. 
Let  R(tu)  =  R(<*0X.  Then  R  is  continuous  at 

R(<u)  c  Ker  L(<d)*J'  =  ^Oq)  +  •  •  •  +RpO0)  • 

This  completes  the  proof. 

Since  is  continuous  at  ujq.  Re  (uj^)  <  0  implies  that  Re  (iu)  <  0 

throughout  some  neighborhood  of  u)q.  Let  (u)o),...,A^  (ujq)  be  the  stable 

1  p 

members  of  a(E (<jjq)  ,A(ujo) ) .  Continuity  at  <jjq  of  and  lemma  6.1  give  the 
following  result. 

Proposition  6.1.  Let  4(uj)  be  the  eigenspace  corresponding  to  the  stable 

P 

eigenvalues  of  the  system  (4.8).  Then  uu-  (ui)  is  continuous  at  uu 

J  j  & 

and  there  exists  a  neighborhood  V  of  such  that  (w)c:A(w)  for 
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every  a>  €V. 

Proposition  6.1  states  that  for  small  perturbations  the  slow 
subsystem  is  at  least  as  stable  when  perturbed  as  it  is  at  'Hq.  Also,  the 
stable  subspace  is  well-behaved  about  iuq. 

We  next  consider  the  controllable  subspace  ft  (uu)  of  the  slow 
subsy s  tem  (4.26). 

Proposition  6.2.  There  exists  R  :  £2"*«(  such  that 

1)  R  is  continuous  at  uj  . 

o 

2)  R(iu)cftg(uj)  VtucQ. 

3)  R(»o).as(«o). 

Proof:  Since 

v(L  (u^B  (u>) )  *  ji  (L  (uj)  )Lv(B  (<u)  ) , 

S  S  a  b 

at  “•  v(Lg  (u))iBg  (iu)  )*  is  continuous  at  ajQ.  Hence,  from  lemma  4.2  there 

exists  R^  :  Q-TTi,  continuous  at  <do,  with  R^Cuj)  c  Im(Lg (<u)1Bg  (cu) )  for  all 

uj  6Q  and  R.  (uu  )  =  Im(L  (ou  )iB  (uj  )).  Let  R  be  as  in  lemma  6.1  with  p  =  r. 
10  s  o  s  o 

Then  1)  holds  and 

R(u»)CR1(ou)  +  ...  +Rr(<u)cImBs(iU)  +In(Ls(u»)Bs(a>))  +  ...  +  Im(Ls  (cu)rBg  (u»)) 

=  ftg  (ou) 

for  all  <u  € Q  so  2)  holds.  3)  follows  similarly.  This  completes  the  proof. 

Thus,  the  slow  subsystem  is  at  least  as  controllable  when  per¬ 
turbed  as  it  is  at  u)  .  The  slow  controllable  subspace  is  well-behaved 

about  <u  . 
o 

To  conclude  this  section  we  consider  stabilizability . 
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Proposition  6.3.  If  (oj  )  is  a  controllable  mode  then  X..  Ciu)  is  controll- 
- 1  i  o  i 

able  throughout  some  neighborhood  of  cu  . 

Proof:  By  hypothesis, 

Im(\ .  (iu  )  I  -  L  (o>  ) )  +  Im  B  (u>  )  =  S  (<u  ) . 

10  so  so  o 

From  proposition  4.1  and  lemma  4.2  there  exist  R^,  :  Q-%  continuous  at 

u>  ,  with 
o 

R1(m)cim  ,(Xi«I  -  Lg(u>)) 

and 

R2(o))clm  v  (Bg(<u)> 
for  all  oi6Q  and  such  that 

W  +  W  3S(lJUo)“ 

From  lemma  6,1  and  the  constant  dimensionality  of  S, 

Im  fit  (\L  (u))I  -  Lg  (<«)  )  +  Im  v(Bg  (id)  )  =  X 
throughout  a  neighborhood  of  Q  and  the  proof  is  complete. 

Corollary:  If  the  slow  subsystem  is  stabilizable  at  tu  then  it  is  stabiliz- 
- 1  '  0 

able  throughout  some  neighborhood  of  uuq. 

6.2.  Fast  Subsystem 

Since  many  different  singularly  perturbed  systems,  with  various 
types  of  fast  mode  behavior,  share  the  same  limiting  descriptor  system,  no 
information  concerning  stability  of  the  perturbed  fast  subsystem  can  be 
extracted  from  the  fast  subsystem  at  the  singular  point.  Hence,  we  cannot 
make  statements  analogous  to  those  of  the  previous  section  about  stabiliz- 
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ability  and  convergence  of  the  stable  eigenspaces.  However,  we  can 
describe  behavior  of  the  controllable  subspace. 

Let  ^(uu)  be  the  controllable  subspace  of  the  fast  subsystem 
(4.27).  If  Lf(ou)  is  not  invertible  at  a  certain  <ju,  then  (4.27)  must  be 
decoupled  according  to  the  descriptor  system  decomposition  outlined  in 
Chapter  2.  The  resulting  slow  and  fast  subsystems  (subsystems  of  (4.27)) 
then  have  well-defined  controllable  subspaces  as  given  in  Chapter  3.  The 
vector  sum  of  the  slow  and  fast  controllable  subspaces  is  then  the  controll¬ 
able  subspace  of  (4.27). 

Proposition  6.4.  There  exists  R  :  ft— 771  such  that 

1)  R  is  continuous  at  uo  . 

2)  R(ou)  Cftf(tu)  V  cu  €  Q  , 

3)  R(wo)  -«f(«u0). 

Proof :  Choosing  an  arbitrary  u)  and  proceeding  according  to  the  algorithm 
(2.10)  -  (2.14)  yields  the  decomposition  F^©F2  =  F(ua)  with 

Lf(u>)  I  Fj_  -  Ll 

and 

L£<»>  1  F2  ■  h 

where  is  invertible  and  L2  is  nilpotent.  In  fact,  the  system  decomposi¬ 
tion  takes  the  form 

x^  =  L^x^  +  L^B^u 


L2x2=x2  +  B2u 


where  B^Pj.  p  Bf(u>)  and  B2  =  Pp^Bf  ((U)  . 


subspace  corresponding  to  subsystem  1  is 


The  part  of  the  controllable 
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1 


=  Im(L^B^)  +  .  .  .  +  Im(L^  B^)  =  (Ira  B^  +  . . .  +  Im  B^) 

=  Im  B^+.«.+  ^Im  B^ 

by  the  Cayley  -  Hamilton  theorem  and  the  -  invariance  of  ft  .  Corresponding 
to  subsystem  2  we  have 

R2  =ImB2+  ...  +L2tl'r'1ImB2. 


Thus , 


flf  (<u)  =  ftL  Sft2 

_  _  „  i 

=  Im  +  ... +L^(id)  ImB^(uu), 


The  result  follows  from  lemmas  4.2  and  6.1.  This  completes  the  proof. 

It  is  gratifying  to  note  that  we  have  defined  controllability 
for  descriptor  systems  in  such  a  way  that  controllability  at  implies 
controllability  throughout  a  neighborhood  of  u>Q. 


6,3,  Application  to  Pole  Placement 

From  propositions  6.2  and  6.4  and  lemma  6.1  it  follows  that 
statements  identical  to  those  in  propositions  6.2  and  6.4  hold  for  il(tu)  = 
(to)  ®R^(uu) .  As  an  application  of  this  and  the  results  of  the  previous 
two  sections  we  now  show  how  some  modes  of  the  perturbed  system  can  be 
placed  approximately  by  designing  a  feedback  gain  for  the  system  at  u)q. 
Here  we  are  generalizing  results  of  [16], 

First,  since  controllability  at  implies  controllability  for 
small  perturbations,  the  existence  of  a  feedback  gain  that  achieves 
arbitrary  eigenvalue  assignment  in  the  perturbed  system  can  in  some  cases 
be  established  by  testing  the  system  at  u)  for  controllability.  If  a  mode 


of  Che  limiting  system  is  controllable  and  if  it  is  shifted  as  desired  by 
linear  feedback  then,  from  local  continuity  of  slow  eigenvalues,  the  same 
feedback  gain  applied  to  the  perturbed  system  results  in  eigenvalues  only 
slightly  different  from  those  desired.  Hence,  for  small  perturbations, 
modes  controllable  at  can  be  approximately  assigned  as  desired  by  consi¬ 
dering  only  the  limiting  descriptor  system. 

In  fact,  as  outlined  in  section  3.5,  given  a  certain  degree  of 

controllability  of  the  limiting  system,  a  feedback  gain  may  be  constructed 

such  that  the  closed-loop  system  at  uoq  has  slow  subsystem  (in  the  descriptor 

sense)  of  dimension  rank  E('JU  )  with  all  modes  controllable.  These  can  be 

o 

assigned  with  linear  feedback  yielding  an  approximate  assignment  in  the 
perturbed  system.  Unfortunately,  this  is  the  best  that  we  can  do.  All 
information  about  the  position  of  the  remaining  n-rank  E(u>o)  eigenvalues  of 
the  perturbed  system  is  lost  at  In  order  to  place  the  remaining  eigen¬ 

values,  the  subsystem  decomposition  must  be  calculated  at  the  perturbed 
value  of  w  and  the  gain  calculated  accordingly.  Nevertheless,  some  computa¬ 
tional  convenience  is  achieved  since  tne  feedback  gains  may  be  calculated 
for  the  slow  and  fast  subsystems  individually.  The  modal  separation 
eliminates  some  of  the  problems  associated  with  stiff  numerical  computations. 
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CHAPTER  7 

THE  LINEAR  QUADRATIC  REGULATOR 


7.1.  Preliminaries 

In  this  chapter  we  consider  the  optimal  control  problem  with 

quadratic  cost  and  singularly  perturbed  system  constraint.  A  similar 

problem  pertaining  to  the  standard  system  (1.5)  was  considered  in  [15]. 

We  shall  need  to  solve  the  regulator  problem  for  the  descriptor  system 

at  The  LQ  regulator  has  been  considered  in  [10]  for  descriptor  systems 

using  dynamic  programming,  but  we  shall  take  a  Hilbert  space  approach.  The 

Hilbert  space  methodology  is  more  suitable  for  dealing  with  questions  about 

convergence  of  the  optimal  control  with  respect  to  uu. 

2 

Let  L  (x)  be  the  set  of  Lebesgue  measurable  maps  x  :  [0,®)  -X 

satisfying 

®  i2 

/  ;|  x(t)  ||  dt  <“.  (7.  1) 

o 

2 

After  identifying  functions  which  are  equal  almost  everywhere,  L  (X)  is  a 
Hilbert  space  with  inner  product  defined  by 

<x,y  >  *  <x(t),y(t)  >dt  (7.2) 

for  all  x,y£L^(X).  Define  L^(U)  similarly.  L^(X)xL^(U)  is  also  a  Hilbert 
space  with  inner  product 


<(x,u)  ,  (y,v)  >- <x,y  >+ <u,v  >  .  (7.3) 

1 

1  Recall  the  definition  (2.21)  of  e(T)  for  T€Hom(X,X).  For 

1  x  6L  (X)  it  is  known  (see  [22], p.  158)  that  if  T  is  stable  then  the  con- 

™  2 

volution  e(T)*x  :  [0,®)  -X  belongs  to  L  (X)  and 

I 

I 


36 


ie(T)*xi|*  !lxijj“  !|  etT||dt  <®. 


(7.4) 


Throughout  this  chapter  we  shall  consider  only  singularly 
perturbed  systems  with  )  =0,  invertible  for  cu^cu^,  and  Lg(u>o) 

stable.  We  further  assume  that  Lf(u>)  is  stable  for  tu^u^.  Hence, 


Tj  ;  [0,®)  -X  defined  by 


tLf(<u) 


-1 


<  CL  <<«) 

e  xo^(cu)+e  ^  xfo(cu)  ^  UJ^(JJ,- 


\<c>  =  \ 


so 


CLs<V 

e  xso(<V 


(7.5) 


if  ou  =  ou 


is  in  l  (X)  for  <jj  in  some  neighborhood  of  cu  .  We  now  restrict  Cl  to  a 

9 

neighborhood  of  ojq  such  that  Tj^^L^X)  for  all  cu  £Q.  Clearly,  Tl^  is  the 
natural  response  of  (4.8). 

If  we  define  a?  :  L^(U)  -L“(X)  by 


e(Ls(o))*  Bg  (u))u  +  e(Lf  (<u)  )*L£(<u)  LBf(uJ)u  if  uj  a>Q 


J  (u)  = 

U) x 


(7.6) 


e(L  (tu  ))*B  (cu  )u  -  BC((JU  )u  if  cu  =  cu  . 
'  s '  o  s  o  f  o  o 


then  ^(u)  is  the  forced  response  of  (4.8). 


7.2.  Problem  Formulation 
Let 

A(cu)  =  l(x,u)  6L2(X)  xL2(U)  !x«a7  (u)i.  (7.7) 

2  2 

Since  is  a  linear  map,  A(<u)  is  a  subspace  of  L  (X)  xL  (U) .  In  fact,  from 

Ju 

(7.4)  it  follows  that  J  is  continuous  so  A(cu)  is  closed  for  all  cu  €Q.  The 
solution  of  the  regulator  problem  minimizes  the  cost  functional 
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J(x,u)  =j  ,x(t)  ;"+  u(t)  !2dt  (7.8) 

o 

2 

over  the  constraint  set  (T.  ,  0)  +  A(ui)  for  each  <JJ  Since  J(x,u)  = '' (x,u) '■  , 

2  2 

we  can  view  the  problem  as  a  minimum  norm  problem  in  L  (X)  xL  (U).  Since 
A(w)  is  closed,  we  have  from  the  projection  theorem  (see  [23])  that  a 
unique  pair  (x^u^)  ^  (~^> 0)  + ''-(UJ)  chat  minimizes  J  exists  for  each  uu 
Furthermore, 

(vV  =  (v°>-pA(uo(V°>  <7"9> 


where  P,  .  is  the  orthogonal  projection  operator  on  A(tu) . 

The  problem  that  this  chapter  addresses  is  not  that  of  finding 
explicitly  the  solution  of  the  regulator  problem  since  this  has  already 
been  done  via  the  algebraic  Riccati  equation.  The  Riccati  equation  will 
not  come  into  consideration  to  any  significant  extent  in  this  chapter.  The 
problem  that  we  shall  consider  can  be  dealt  with  much  more  simply  from  a 
geometric  point  of  view. 

Our  problem  is  that  o'  establishing  conditions  under  which 
2 

x  —x  and  u  -*u  in  the  L  sense  as  uj—io  .  If  (x  ,u  )  converges  then  for 
a;  u>  a>  uu  o  oi  yu 

o  o 

each  S>0  there  exists  a  neighborhood  V  of  id  such  that  if  u>  £v  then 

o 


and 


(7.10) 

(7.11) 


o  o 


(7.12) 


since  and  J  are  continuous.  Thus,  for  small  perturbations  about  uo  , 

the  solution  of  the  regulator  problem  at  uu  can  be  applied  to  the 
o 
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perturbed  system  with  only  a  slight  deviation  from  the  optimal  trajectory 

and  optimal  cost.  It  will  be  shown  that  solving  the  problem  at  tu 

explicitly  for  is  simpler  computationally  than  solving  the  problem  at 
o 

for  u^.  Hence,  a  computational  savings  can  be  achieved  at  the  cost 
of  slight  suboptimality. 

7.3.  Reduced -Order  Solution  at  u) 

_ _ _ o_ 

As  a  result  of  the  reduction  in  order  that  occurs  in  the  system 

at  a  .  solving  the  regulator  problem  explicitly  at  cu  turns  out  to  be 
°  o 

simpler  than  solving  it  for  some  other  U).  To  see  tnis,  choose  an  ortho¬ 
normal  basis  (eL,...,er)  of  S(ouq),  an  arbitrary  basis  (er+^>  •  •  •  >en)  of 


F<V 

and  an  orthogonal  basis  (v^,... 

’  Vm 

)  of  U  with  llvjl  => 

■Y,  i  —  l,...,m 

for  some  Y  >0.  If 

u  €  L2  (U  )  and  x  €  L2 

(X) 

with  x  *  x  +  x-,  : 

s  f  ’ 

Ks(t)  es(u)o). 

xf(t) 

€F(mo),  and 

xs(t) 

=  + 

•  *  * 

+  “r(t)er 

(7. 

.13) 

Xf  (t) 

- 

+  . 

. .  +  8(  t )  e 
n-r  n 

(7. 

.14) 

U(t) 

=  Yl(C)Vl+ 

+  Y  (C)v 
mv  '  m 

(7. 

,13) 

then 

J(x,u)  -J 

* 

a(t) 

o 

A 

a(t)  +  2ar(t) 

n3( 

:o  +  @(t)*QP(t)+Y 

2Y (t)*Y (t)dt  (7. 

16) 

where  af(t),  8(t),  and  Y(t)  are  column  vectors  consisting  of  the  c^Ct), 
8^(t),  and  Y^(t)  and 
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N 


Q  = 


<e,,e  >  ...  <  e,  ,e  > 

1  r+1  1  a 


<  e  ,e  ,>...  <e,e  > 
r*  r+1  r*  a 


<  e  , , ,e  .  >  . . .  <e  ..  ,e  > 
r+1’  r+1  r+1*  n 


<e  ,e  >  ...  <e,e> 

ir  r+1  n*  n 


(7.17) 


(7.18) 


If  (x,u)  e  (Tl^.O)  +A((Uo)  then  x(C)  =xg(t) +xf  (t)  =xs(t)  -  Bf(0UQ)u 
so  we  may  equivalently  minimize 


J(xs,u)  =  J(xs  -  Bf  (ujo)u3u)  »j  a(t)  a(t)  -  2a(t)  NKY(t) 

+  Y  ( t  )*  (Y2I+K*QK)f(t)dt  (7.19) 

where  K=MatB,(<u  ). 

f  o 

In  [31],  pp.  46-48  it  is  shown  that  the  optimization  problem 
with  cost  (7.19)  and  system 

or(t)  -Gcr(t) +HY(t)  (7.20) 

may  be  reduced  to  that,  of  minimizing 

j(a,0)  =/  0f(t)*(I  -NK(Y2I  +  K* QK)_1TC* N*  )a(t)  +  0(t)*  (Y2I  +  K  QK)0(t)dt  (7.21) 


subject  to 


d(t)  =  (G  -H(Y2I  +K* QK)-1K* N*)a(t)  + 


H0(t) 


(7.22) 


if  I-NK(Y2I+K  QK)  LK  N  is  positive  semidefinite. 

o  ^  —  I  ★  Jc 

Lemma  7.1.  There  exists  Y  >0  such  that  I  -  NK(Y  I+K  QK)  K  N  is  positive 
definite. 

Proof:  The  matrix  T»NK(Y  I+K  QK)  K  N  is  clearly  Hermitian,  positive 


semidefinite  since  Q  is  positive  definite.  We  need  to  show  only  that  for 
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some  Y  >0  Che  eigenvalues  of  T  are  less  than  unity.  Since  the  eigenvalues 
of  a  matrix  are  bounded  above  by  its  norm,  it  is  sufficient  to  show  that 

:<y2l+K*QK)'li;,<^ 

for  some  Y . 

According  to  a  well-known  result,  there  exists  p>0  such  that 

!Il!I  <  p  5(l) 

for  all  LcC™™  where  i  is  the  spectral  norm.  Let  r*  be  the  minimum  eigen- 
* 

value  of  K  QK  and. let 

Y  >VmaxCp!iNKi|2  oJ . 

Then 

ii  (Y2x  -hk'qk) -1j{  <P  i((Y2I  +  K*QK)'1>  -  ~  . 

Y^  +  M-  ,'lNKir 

This  completes  the  proof. 

Setting  G=Mat  Ls(yaQ)  and  H=Mat  Bs  ((jjq)  and  choosing  Y  as  in 
lemma  7.1  we  may  solve  the  regulator  problem  at  by  minimizing  J  subject 
to  (7.22)  which  necessitates  the  solution  of  an  rth  order  Riccati  equation. 
Compared  to  the  nth  order  problem  for  the  rth  order  problem  is  a 

considerable  computational  simplification.  The  remainder  of  this  chapter 
is  demoted  to  finding  conditions  under  which  the  solution  of  the  reduced 
order  problem  is  close  to  that  of  tne  full  order  perturbed  problem. 

7.4.  Convergence  of  A(ua) 

In  order  to  establish  convergence  of  the  optimal  solution  (x^uj 
it  is  convenient  to  establish  convergence  of  A (oo)  with  respect  to  a  certain 
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I 

I 

I 

I 


2  2 

metric  on  the  closed  subspaces  of  L  (X)  xL  (U).  Let  ‘ff(  be  the  set  of  closed 
subspaces  of  a  Hilbert  space  H.  In  [32]  and  [33]  a  metric  P  on  /)(  is 
discussed  giving  two  equivalent  expressions.  Let 

sup  sup 

p(R,T)  -  ,|P  -  P  =maxl  x€R  D[x,T],x€T  D[x,R]i  (7.23) 

llxll  s  1  ilxii  £  1 

for  any  R,  T  657;  where 


D[x,T]  -JgiU-yl! 


(7.24) 


Note  that  (/?!, p)  is  a  generalization  of  the  metric  space  considered  in 


Chapter  4. 


-1 


•°°  i  tLf 

Theorem  7.1.  If  J  lie  ||dt  -*0  as  cju-»(jo  then  the  mapping  u)-A(uj)  is 

o 

continuous  at  with  respect  to  P. 


Proof;  observe  that 


sup  D[  (x,u),A(uu  )]2  £  sup  D[(»<  (u),u),A(oj  )]2 
(x,u)eA(ou)  °  u(L^(U) 


|i(x,u)||  s  1 


u  S  1 


sup  0  inf  ,  iK^(u)  (u),u-v)||2  ssup  2  I'k^u)  (u);|2 


u  €L2(U)  v  £L2(U) 

HI s  1 


u  6  L2  (U  )  ^ 

'lull  si 


j  \o/  -  g/ 

,rW  U) 


Simi larly. 


sup  D[(x,u),A(uj)]  s  ))✓  -J 


(x,u)  feA((u  ) 
il(s,u)||Sl° 

®  tLf(uu) 

so  from  (7.4) , convergence  of  J  i|e 

o 

guarantee  that 


uu  yj„ 


-1 


i|dt  to  0  and  stability  of  Ls(tuo) 
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P<A(*),A(<do))  -0 


as  UJ_’UJ0  and  the  proof  is  complete. 


7.5.  Convergence  of  the  Optimal  Solution 

We  now  give  sufficient  conditions  under  which  tu-(x  ,u  )  is 

'  UT  Ur 

continuous  at  a  ,  , 

tLf^  ,,  ;  1 2 

Theorem  7.2.  If  J  lie  ildt-OandJ  j|e  ,!  dt-*0as  io ~ ' to  then 


a>  —  (x  ,u  )  is  continuous  at  . 
tu  0)  o 

Proof :  First  observe  that 


•I. 


UO  (JU 


-i  jz 

TL  llAfr  He  f  ||2! 


|xfo(oj)  H  dt  +  l|e(Ls(aj))xso(m)  -  e(Lg  (U,0>)XS0(U>0)  I 


so  stability  of  L  (a;  )  and  tne  second  hypothesis  of  the  theorem  imply 
s  o 

-  But  the  first  hypothesis  guarantees  that 
o 

PA(u))  "*  PA(uJq) 

from  (7.23).  From  (7.9)  we  have 

'  (xoi  ’um  )i'S’llVT|<u  ^  +  l'K(<v) 'PA(w  )^T|a>  '' 

0  0  O  0  0 

and  the  desired  result  follows. 


Since  the  hypothesis  of  theorem  7.2  implies  that 

r“l  =i-£<u,)-1 

J  lie  II  dt  <» 

o 


(7.25) 


throughout  a  neighborhood  of  a)  and  since 

1  0 

=  tL  (a))"1  2tRe  FTij"  -  if  Re 

J  ' |e  £  il  dt  ij  e  J  dt  =  < 
o  o 


g.  (u>) 


a  0 


v  2  Re 


aj(a») 


if  Re 


(7.26) 


a.  («> 


<  0 
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for  all  uj  Sfl,  j«l, . . .  ,n-r,  it  follows  that  it  is  necessary  for  L^(uj)  to  be 
stable  throughout  a  neighborhood  of  in  order  that  the  sufficiency  condi¬ 
tion  of  theorem  7.2  hold. 

The  sufficient  condition  of  theorem  7.2  guarantees  that  the 

reduced  order  problem  at  <Do  may  be  solved  yielding  only  a  slightly  sub- 

optimal  control  which  is  close  to  the  optimal  control  and  which  generates  a 

2 

trajectory  close  to  the  optimal  trajectory  in  the  L  sense.  Under  most 
circumstances  it  is  reasonable  to  interpret  the  convergence  criteria  of 
theorems  7.1  and  7.2  as  conditions  on  the  behavior  of  the  boundary  layer 
in  singularly  perturbed  systems,  for  if  the  natural  response  converges 
uniformly  to  zero  in  [€,<*>)  the  integrals  are  essentially  measures  of  the 
intensity  of  the  boundary  layer  effect  in  [0,€).  The  integral  convergence 
conditions  state  that  the  effect  of  the  boundary  layer  becomes 
vanishingly  small  as 
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CHAPTER  8 

ALTERNATIVE  FORMULATIONS  AND  CONCLUSIONS 

8.1.  Algebraic  Interpretations 

Although  it  was  not  explicitly  stated  in  previous  chapters,  there 
are  abstract  algebraic  interpretations  for  many  of  the  results  encountered 
so  far  in  our  study  of  singularly  perturbed  systems.  These  interpretations 
have  not  been  explored  in  much  detail  yet,  but  they  are  presented  here  for 
completeness  and  as  a  suggestion  for  further  study. 

To  begin  with,  consider  the  set  of  all  mappings  \  :  Cl  — C,  continuous 
at  Such  mappings  have  been  considered  extensively,  starting  with  Chapter 

4,  but  we  have  not  considered  properties  of  the  set  of  all  such  maps.  Denote 
the  set  by  C(wq).  With  little  effort  it  can  be  shown  that  C(ioq)  is  a  commuta¬ 
tive  ring  with  identity  using  pointwise  addition  and  multiplication.  Letting 

2(.u  )  denote  the  subset  of  C(oj  )  consisting  of  all  X  with  \(w  )=0,  it  can  be 
r  o  o  o 

shown  that  is  an  ideal  of  C(Wq).  Lemma  4.3  may  be  interpreted  as  a 

factorization  theorem  for  polynomials  over  C(u)q). 

Let  ax  be  the  set  of  maps  x  :  Q  — X,  continuous  at  Using 

pointwise  addition  and  scalar  multiplication,  Qx  is  a  C(u)q)  -  module.  Let 
7H,  be  the  set  of  subspaces  of  X  and  R  :  Q-#i  be  continuous  at  with  dimen¬ 
sion  p  for  all  co  SCI.  For  example,  R  may  be  the  slow  or  fast  subspace  map 
S  or  F  as  defined  in  Chapter  4.  R  may  be  identified  with  a  submodule  of 
Gx  in  the  following  natural  way.  As  outlined  in  lemma  4.2  choose  a  basis 
(x^(rn) , . . .  ,Xp(m))  of  R(co)  with  x^  continuous  at  Then  x^  £Qx  and 

lXp...,Xpi  is  linearly  independent  since 
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a1(uJ)x1(u>)  +  . . .  +ap(u>)xp(cu)  =  0 

(8.1) 

implies  crCui)  »0  for  all  u),  i*l,...,p.  Let  (y^Ou) , . . .  ,yp 
basis  of  R(u))»  continuous  at  <uq.  Then 

(u>))  be  another 

y1(<«)  *  P1(w)x1(u))  +  . . .  +  @p0u)xp(u0 

(8.2) 

with  9  ^C((juo>  by  lemma  4.1.  Hence, 

spanly^  . . .  ,ypJ  c  spanix^  . . .  ,xp} . 

(8.3) 

By  reversing  the  argument, 

span[x^, . . . ,xpj c  spaniy^, . . . ,ypj  (8.4) 

so  we  may  naturally  and  without  ambiguity  identify  R  with  spanlx^, . . . ,xpj . 

In  our  study  of  singularly  perturbed  systems  we  considered 

operator  valued  maps  A€HR(X).  Members  of  HR(X)  can  be  identified  with 

linear  transformations  on  the  submodule  R  by  setting 

(Ax)(ui)  *  A(uu)x(uj). 

(8.5) 

H  (X)  admits  the  structure  of  a  C(cu  )  -  algebra.  Let 

K  0 

(A  +  B)(U)  ■  A(W)  +BIW) 

(8.6) 

(AB)(<U)  «  A(<u)B(U» 

(8.7) 

(U)  «H)  -  \<«i)A(<u) 

(8.8) 

for  A,  B  fcl^(X). 

Following  Che  same  line  of  reasoning  ic  can  be  seen  Chat  Che  sec 

<1.  of  all  maps  u  :Q“*U  Chac  are  conCinuous  ac  a)  is  a  C(u>  )  -module.  Also, 
u  O  O 

Hr(U)  is  a  C(ujo)  -  module  of  linear  transfonnacions  from  <2^  inCo  Q^. 

Consider  pencils  (i.e.  first  degree  polynomials)  over  the  algebra 


H  (X).  Using  lemmas  4.1  and  4.2  it  is  easily  shown  that  all  the  bases  of 

the  submodule  R  have  the  same  number  of  elements.  For  G,  L£H^(X)  each 

basis  of  R  determines  rxr  matrix  representations  of  G  and  L  with  entries 

in  C(u)  ).  It  also  determines  a  matrix  representation  of  the  pencil  (G,L)^ 
o 

with  entries  in  the  C(<u  )  -  algebra  C(tu  )[s]  of  polynomials  over  C(U)  )  in 

o  u  U 

the  indeterminate  s.  The  determinant  of  the  pencil  (E,A)  may  be  defined 
by  forming  the  determinant  of  Mat  (E,A)  in  the  usual  way  with  respect  to 
some  basis  yielding 

det(E,A)  €C(U>  )[s] .  (8.9) 

A  simple  argument  shows  that  det(E,A)  is  independent  of  the  basis  chosen. 

Consider  A  £H  (X)  and  let 
K 

p 

det(I,A)  -Y  i21  (s-i^)  (8.10) 

where  I  is  the  identity  element  of  H  (X)  and  Y  £C(uj  )  is  invertible. 

(Y  is  invertible  if  and  only  if  Y(w)?t  0  for  all  a  €  Q) .  Define 

a(A)  =  iT(ili=l,...,p}  (8.11) 

The  7|^  can  be  considered  as  eigenvalues  of  A. 

We  are  now  in  a  position  to  interpret  the  central  singular  pertur¬ 
bation  decomposition  result,  theorem  4.1,  algebraically.  Suppose  that 
E,  A(Hx(X)  with 

r  n-r 

det(E,A)  -0o(1"1(3-V)(i^1(V-l))  (8.12) 

where  6C(ojo)  is  invertible  and  onS^uu^).  This  can  be  done  according  to 

^The  pencil  (G,H)  is  often  written  Gs  -H. 


97 


lemma  4.3.  Theorem  4.1  states  that 

(E,A)  =  N  *  (K,L)  (8.13) 

where  N  is  an  invertible  constant  polynomial  over  H  (X) ,  denotes 

A 

polynomial  multiplication,  and  K  and  L  are  both  S-  and  F-invariant  with 


K  1  S=I  6H  (X) 

J 

(8.14) 

L  1  F  *  I  6Hf(X) 

(8.15) 

c(kIf)  = [a J i=l, . . . ,n-r) 

(8.16) 

ct(lIs)  = lx  li»l,...,r}. 

(8.17) 

Theorem  4.1  may  be  interpreted  as  a  canonical  factorization  result  for 
regular  pencils  (i.e.  with  det(E,A)  ^  0)  over  H^(X). 

8.2.  Geometry  of  the  Space  of  Linear  Systems 

Let  T  be  the  complex  Euclidean  space  of  ordered  pairs  of  nxn 
matrices  (E,A).  There  is  an  obvious  one-to-one  correspondence  between  F 
and  differential  equations  Ex=Ax.  Hence  each  point  of  F  can  be  inter¬ 
preted  as  a  linear  system  of  one  of  the  following  three  types:  1)  a  state 

variable  system  if  E  is  nonsingular,  2)  a  descriptor  variable  system  if  E 

? 

is  singular  and  det(Es-A)  $0,  3)  a  degenerate  system  if  det(Es-A)  a0. 

Viewing  linear  systems  in  this  way  is  natural  since  a  small  perturbation  of 
a  given  system  in  the  Euclidean  norm  is  equivalent  to  a  small  perturbation 
in  the  system  parameters. 

Studying  the  geometry  of  F  adds  valuable  insight  into  the  nature 
of  descriptor  variable  and  singularly  perturbed  systems.  With  little  effort 
| 

This  implies,  of  course,  that  E  is  singular. 
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it  can  be  shown  that  descriptor  and  degenerate  systems  together  form  a 
hypersurface  in  F  contained  in  the  boundary  of  the  set  of  state  variable 
systems.  The  property  of  being  state  variable  is  a  generic  property  of  F. 

It  is  an  unfortunate  fact  that  for  every  descriptor  system  in  F 
there  exists  a  sequence  of  state  variable  systems  converging  to  it  with 
the  corresponding  sequences  of  eigenvalues  diverging  to  +00  .  That  is,  if 
a  descriptor  system  is  perturbed  in  the  wrong  direction,  the  peiturbed 
system  will  have  tremendous  instability,  the  smaller  the  perturbation,  the 
greater  the  instability.  It  is  as  if  every  descriptor  system  is  perched 
precariously  on  the  edge  of  a  cliff.  A  step  in  one  direction  will  result 
in  only  a  slight  change  in  its  characteristics.  A  step  in  the  other  direc¬ 
tion  will  have  disastrous  consequence.  The  importance  of  establishing 
simple  conditions  that  guarantee  trajectory  convergence  is  clear.  If  a 
designer  fails  to  account  for  the  possibility  that  a  descriptor  system's 
parameters  are  slightly  different  from  what  he  thinks  they  are,  his  whole 
design  could  fail  miserably. 

Let  DCT  be  the  set  of  degenerate  systems.  One  way  to  view  the 
question  of  trajectory  convergence  is  to  consider  the  map  Sx  :  r-D—£o(x) 
which  associates  with  each  pair  (E,A)  the  solution  of  Ez=Az  for  initial 
condition  x  €X.  If  a  topology  is  placed  on  ^Q(X)  (e.g.  see  Chapter  5)  then 

we  need  to  ask  questions  about  the  weakest  topology  on  T  that  makes  §x  continuous . 

Once  the  nature  of  the  resulting  neighborhoods  of  a  descriptor  system  is  understood, 

the  behavior  of  the  solutions  of  a  singularly  perturbed  system 

(which  is  nothing  more  than  a  map  from  C2  into  T,  continuous  at  u>  or  the 

parameterization  of  a  particular  path  in  F)  can  be  determined  by  checking 

to  see  if  an  arbitrarily  small  neighborhood  of  the  descriptor  system 
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contains  the  image  of  some  neighborhood  of  . 

So  far  in  this  section  we  have  considered  only  pairs  of 
matricies  corresponding  to  unforced  systems  with  fixed  initial  conditions. 

In  order  to  study  the  behavior  of  a  forced  system  with  a  parametrically 
varying  initial  condition,  T  must  be  the  space  of  4-tuples  (E,A,B,x)  where 
B  is  nxm  and  x  is  n  x  1.  Our  previous  discussion  of  geometry  and  induced 
topology  carries  through  with  only  minor  changes. 

8.3.  Suggestions  for  Further  Research 

In  the  area  of  descriptor  variable  theory  there  are  many  avenues 
which  have  yet  to  be  explored.  For  example,  although  observability  of  descriptor 
systems  has  been  considered  in  [7],  the  descriptor  variable  equivalent  of 
observers  from  state  variable  theory  have  not  been  developed.  In  a  stochastic 
environment  the  Kalman  filter  might  also  have  a  natural  extension  to  descriptor 
systems.  There  are  many  fundamental  control-theoretic  concepts  £  ich  as  the 
Maximum  Principle  and  various  stochastic  and  adaptive  control  techniques  that 
have  yet  to  be  considered  in  the  context  of  descriptor  variable  theory. 

Clearly,  the  problem  of  trajectory  convergence  in  singularly 

perturbed  systems  has,  for  tne  most  part,  not  been  solved  except  for  certain 
standard  systems.  The  general  case  that  we  have  considered  still  requires 
a  great  deal  more  work.  The  resolution  of  this  issue  is  essential.  As  we 
have  seen  in  the  last  section,  the  survival  of  the  descriptor  variable 
approach  to  system  modeling  depends  on  it. 

As  with  descriptor  systems,  many  system-theoretic  concepts  have 
not  as  yet  been  extended  to  singularly  perturbed  systems.  In  Chapter  7  we 
studied  the  regulator  problem  for  Ls (<Uq)  stable  and  L^u)  )  =0.  The  same 
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problem  needs  to  be  considered  with  the  two  assumptions  dropped.  It  is 
safe  to  say  that  if  a  problem  has  not  been  studied  in  the  context  of 
descriptor  systems  then  it  needs  to  be  studied  in  the  context  of  generalized 
singularly  perturbed  systems.  If  a  given  system  is  close  to  a  descriptor 
system  in  the  Euclidean  norm  then  conditions  are  needed  to  insure  that  one 
need  only  consider  the  nearby  descriptor  system.  If  a  designer  is  guaranteed 
that  the  application  of  some  design  technique  to  the  descriptor  system  will 
yield  results  close  to  those  that  would  come  from  working  with  the  given 
system,  then  he  may  choose  to  apply  that  technique  to  the  reduced  order 
descriptor  system.  Such  an  action  often  results  in  increased  computational 
efficiency.  Of  course,  the  price  is  always  inferior  system  performance. 

8.4.  Conclusions 

In  this  thesis  three  central  points  have  become  clear.  First, 
there  are  many  alternative  ways  to  view  descriptor  variable  and  singularly 
perturbed  systems.  They  range  from  the  matrix  oriented  approaches  which 
exist  in  most  of  the  literature  to  the  geometric  theory  developed  in  Chapters 
1  through  7  to  the  algebraic  ideas  discussed  briefly  in  this  chapter. 
Certainly,  there  are  other  interpretations  as  well  that  no  one  has  even 
thought  of  yet.  The  more  ways  that  exist  to  look  at  a  problem,  the  more 
likely  it  is  that  the  problem  will  be  solved  in  the  near  future. 

The  second  point  is  that  descriptor  systems  must  be  considered 
as  members  of  the  space  F.  Since  they  are  located  in  such  precarious 
positions  in  F,  failure  to  consider  their  spatial  relationships  with  nearly 
state  variable  systems  could  result  in  unexpected  system  behavior,  to  put 
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it  mildly.  The  question  of  trajectory  convergence  has  yet  to  be  answered 
satisfactorily. 

Finally,  there  are  still  many  important  control-theoretic  con¬ 
cepts  that  have  not  been  extended  to  singularly  perturbed  systems.  We 
have  made  some  progress  in  the  pole  placement  and  regulator  problems  in 
Chapters  3  and  7,  but  many  other  problems  still  exist. 


s*v  v» 
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